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ABSTRACT 

We present new mass measurements for the supermassive black holes (SMBHs) in the centres of 
three early-type galaxies. The gas pressure in the surrounding, hot interstellar medium (ISM) is 
measured through spatially resolved spectroscopy with the Chandra X-ray observatory, allowing the 
SMBH mass (Mbh) to be inferred directly under the hydrostatic approximation. This technique 
does not require calibration against other SMBH measurement methods and its accuracy depends 
only on the ISM being close to hydrostatic, which is supported by the smooth X-ray isophotes of 
the galaxies. Combined with results from our recent study of the elliptical galaxy NGC 4649, this 
brings to four the number of galaxies with SMBHs measured in this way. Of these, three already 
have mass determinations from the kinematics of either the stars or a central gas disc, and hence join 
only a handful of galaxies with Mbh measured by more than one technique. We find good agreement 
between the different methods, providing support for the assumptions implicit in both the hydrostatic 
and the dynamical models. The stellar mass-to-light ratios for each galaxy inferred by our tech nique 
are in agreement with the predictions of stellar population synthesis models assuming a Kroupa initial 
mass function (IMF). This concurrence implies that no more than '^10-20% of the ISM pressure is 
nonthermal, unless there is a conspiracy between the shape of the IMF and nonthermal pressure. 
Finally, we compute Bondi accretion rates (Mbondi), finding that the two galaxies with the highest 
Mbondi exhibit little evidence of X-ray cavities, suggesting that the correlation with the AGN jet power 
takes time to be established. 

Subject headings: Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: ISM — galaxies: 
individual (NGC1332, NGC4261, NGC4472, NGC4649)— black hole physics 



1. INTRODUCTION 

Supermassive black holes (SMBHs) with masses rang- 
ing from less than a million to a few billion times the 
massi of the Sun, believed to be ubiquitoua in galaxy 



fc Ferrarese 2001; Tremaine et al. 2002; Ferrarese & Ford 



2005 ), and a number of authors have discussed possible 
deviations from the universal, log-linear form typ i cally 



bulg es , arc increasingly being recognizins: as essential in -- 
gredi e nts in the formation o f galaxies (e.g. ^ilk fc Rees 
1998| ; |Di Matteo et al.||2005| ; [Hopkins ct al||2006D . CmT 
pelling evidence for the close link between the formation 
of SMBHs and stellar spheroids is provided by the tight 
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could have an important effect (Giiltekin ct al 

Despite their central importance for understanding 
galaxy evolution, only ~40-50 precise, reliable SMBH 
measurements span '-^S orders of magnitude in Mrh ( Fer 



central velocity dis 


persion, cr* (Gebhardt et al. |2000a 


Ferrarese & Merritt 


20001; Kormendy fc Richstone 1995 


Marconi fc Hunt 2003; Graham et al, |2001). Since ditfer- 



lead to quantitatively different behaviour, both the shape 
of these relations, particularly that between Mbh and cr*, 
and the scatter about them have become crucial diag- 
nosti cs for cons t raining our picture of galaxy forniation 



rarese &: Ford||2005| ; |Graham| pOO^ ; piiltekin ct aL||200S| y. 

These are all derived from the kinematics of stars or a ro- 
tating central gas disc in the deep potential well near to 
the black hole. In external galaxies, kinematical informa- 
tion is mostly restricted to the distribution of velocities 
along the line of sight, so the use of either tracer implies 
understanding the three-dimensional velocity field. Gas 
disc modeling generally assumes a circular velocity field, 
so that deviations from Keplerian motion can give rise 



to errors in Mbh (e.g. [Cappellari et al.|[200^ ). Stellar 



(e.g. King 2003; 


Nipoti et al. 


20031; iGranato et al. 


2004; 


Robertson et al. 


2006). Unfortunately, controversies per- 



kinematics are complicated by the degeneracy between 
the enclosed mass and the stellar orbital structure. To 
overcome this degeneracy, it is standard to use the full 
velocity profile, by either including higher order ni oments 
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( van der Marel fc Franx 1993 ; Bender et al 
ing non-parametric fitting ( [Ccbhardt et al 



19941 ) or us- 
pOOObl ). More 

importantly though, it is essential to minimize assump- 
tions about the stellar orbital structure. For an axisym- 
metr ic system, this genera l ly requires orbit-based models 
fe.g. Gebhardt et al, 2003; van der Marel et alj|l998|; Ivan 
H 



den Bosch et al. SUtiS ). However, these mode 



s are not 
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free of systematic conce r ns, such as whether to regu lar- 
ize ( IVanuri et "Z] |2004|: Jvan der Marel et al.| |1998| ) or 



space coverage (Vahuri et al 



not (e.g. Cebhardt ct al |2003 ) , ens uring adequate 



lase 

j004) and the possible im - 
plications of triaxiality (e.g. van den Bosch et al. 2008). 
Since the inclination of an elliptical galaxy is generally 
poorly known this can introduce significant uncertain- 
ties into the assumed three-dimensional light distribu- 
tion, and cons equently the inferred gravitating mass (e.g. 

avazzi 2005). A reliable recovery of Mbh also depends 
on the accurate modelling of the total gravitating mass 
profile over a wide radial range, so that the results can 
be sensitive to the treatment of the dark matter com- 
pone nt. Unfortunately, most published SMBH measure- 



ment s from stellar dynamics omit this mass component, 
leading to a syst ematic underestimate of Mr h by as much 
as a factor -2 ( pebhardt fc Thomas] p009| ). 

Clearly further progress demands independent verifica- 
tion of the different mass determination techniques. An 



effec ive way to address this question is to compare mul- 
tiple Mbh measurements for the same SMBH made in 
different ways. However, reliable measurements of the 
same black hole by both dynamical methods are rare 
and comparisons between them yield mixed results. To 



date t he most detailed exa.mples are for NGC 4258 ( siopis 
et al| |200£ ; Miyoshi et al. 1995 ), for which the measure- 



ments were marginally discrepant (at 2.6a, but within 
~15% of each other), Cen A, for which both techniques 
agree within a factor ~4 but with systematic discrepan- 
cies between different measurements (even those made 



tistical errors f 


Silge et al.l 


2005: 


Niarconi et al. 


12006; 


Neu- 


maytr et al. 


2007 


: jCappel 


ari et al. 


2009 


), and NGC 


3379, 



for whic h there is general agreement for axisymm etric 
models ( [Shapiro et aL||2006|; [Gebhardt et"aI||2000bD, but 
not with triaxial models ( [van den Bosch et al. 2009) A 
handful of similar comparisons are less compelling due 
to deviations from Keple rian gas motion, or overly sim- 



ple s tellar modelli ng (e.g. Cappellari et al, 
Kleij|n et alj ^002|) 



2002 
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A new way of detecting SMBHs , appropriate for early- 
t ype galaxies, was proposed by Brighenti fc Mathews 



(1999), who pointed out that the gravitational influence 
of a quiescent SMBH on the hot. X-ray emitting inter- 
stellar medium (ISM) should affect the temperature and 
density of the gas sufficiently to be detectable in Chandra 
observations of nearby galaxies. Although their models 
predict the ISM to be globally inflowing, the inflow is ex- 
pected to be highly subsonic so that the ISM should re- 
main very close to hydrostatic equilibrium, enabling Mbh 
to be measured directly with the stand ard hydrostatic X- 



ray m a ss determination meth ods (e.g. GastaldcUo et al 
^007b|; [Humphrey et al.||2006|, hereafterir06[) . If the ii^ 



flow is adiabatic over scales which can be resolved with 
Chandra, the resulting gas compression should produce a 
measurable central temperature spike near to the black 
hole. Unlike other techniques that do not rely on dy- 
namical modelling, such as using the variability or lu- 



mino s ity of an active gala c tic nucleus ( Markowitz et al . 
^003[ ; [O'NeiU et al[ [2005[ ; [Ferrarese fc Forc^ [2005[ ,""Mld 



references therein), hydrostatic Mbh measurements do 
not require calibration against other methods and their 
intrinsic accuracy on a case-by-case basis is limited the- 
oretically only by the extent to which the hydrostatic 



approximation holds (in addition to further simplifying 
assumptions usually adopted to make the problem more 
tractable, and provided there is sufficient spatial resolu- 
tion to resolve the gas properties close to the SMBH). 

Since the local freefall timescales in the centres of 
early-type galaxies are generally much shorter than the 
gas cooling time, hydrostatic equilibrium is expected to 
be quickly established unless the ISM is being actively 
stirr ed up, for example by AGN activity or galaxy merg- 
ing ( [Mathews fc Brighenti 2003). Nonthermal pressure 
due to gas turbulence is not expected to perturb this 
equilibrium substantially; hydrodynamical simulations of 
structure formation suggest i t contributes no mo r e than 



25% of th e total pressure (Evrard et al. 199G : Nagai 



et al. 2007) and possibly significantly less ( l^'ang et al 
2009). This is consistent with the general agreement be- 
tween the masses of morphologically relaxed galaxy clus- 
ters meas ured by X-ray and gravitational le ns ing meth- 



et al 



ods ( e .g. [AUenl [1998[ ; 



Mahdavi et al. 2008). Churazov 



(2008) compared stellar dynamics and X-ray mass 



measurements in two early-type galaxies with moderately 
disturbed X-ray isophotes and, similarly, found that non- 
thermal support contributed no more than ~10-20% of 
the global pressure. The excellent agreement between 
the distribution of dark matter in galaxies, groups and 
clusters inferred from X-ray studies and that predicted 
by cosmological simulations provides f urther corrobora- 
tion of the hydrostatic approximation ( Lewis et "al] [200S ; 
Vikhlinin et al] ^006[ ; [H06[ ; [Buote et al] ^007[) . Never- 
theless, the important question of whether the gas in the 
galaxies studied in this paper is c lose to hydrostatic is one 
we will discuss in detail in § 8.3. 

In a recen t Chandra study oFthe g iant ellipti cal ga laxy 
NGC 4649 ([Humphrey et al] [2008| , hereafter ^408[) , we 
used this technique for the first time to constrain the 
mass of the SMBH. This galaxy has an extremely round, 
relaxed X-ray morphology, indicating little or no large- 
scale gas disturbance, and suggesting little turbulence 
that would perturb hydrostatic equilibrium by no more 
than a few percent (Brighenti et al. 2009| ). Using the 
X-ray data, we constrained Mbh to ~ 3 x 1O^M0, which 
is close to the me asurement from stellar dynamics ( [Geb 



hardt et al. 2003). This agreement provides compelling 
evidence not only that the gas must be close to hydro- 
static, but also that the stellar modelling of this system is 
accurate, effectively ruling out pathological orbital struc- 
ture. 

In this paper, we extend our study to include three 
more objects, bringing to a total of four the number of 
early-type galaxies in which Mbh has been determined 
from hydrostatic X-ray methods. Of these, three already 
have Mbh determined by stellar dynamics or gas kine- 
matics while, for the remaining galaxy, we are presenting 
the first direct Mbh measurement. In § || we discuss the 
sample selection, in § ^ we describe the X-ray data re- 
duction, in § ^ we discuss the deprojection of the stellar 
light, an essential step in modelling the mass, in § |^ we 
outline how we infer the mass from the X-ray data and in 
§ ^ we report the best-fitting results. We discuss possi- 
ble sources of systematic uncertainty in § ^ and reach our 
conclusions in § |8[ All errors quoted in the text of this 
paper, unless otherwise stated, correspond to the 90% 
confidence regions (which, for our Bayesian analysis, im- 
plies the region of marginalized parameter space within 



3 



which the integrated probabihty is 90%). In the figures, 
we typically use 1-a errors (68% enclosed probability), 
for clarity. 

2. TARGET SELECTION 

To select our targets for this study, we searched the Hy- 
perLEDA database for early-type galaxies (T<-1) with 
heliocentric recessional velocities < 2300km s~^ (corre- 
sponding to a distance '^33 Mpc) and central velocity 
dispersion, ct* > 280km s~^, (corresponding to a black 
hole mass 5 x lO^Mo; Tremaine et al. 2002). Restrict- 
ing ourselves to gas-rich galaxies that have been observed 
with Chandra^ our sample contained 9 objects. We omit- 
ted three galaxies (M 87, NGC 3998 and IC 1459) as they 
contain X-ray emitting AGNs which may make ISM 
temperature measurements at < 2" scales challenging. 
This left 6 objects: NGC 1332, NGC 1399, NGC 4261, 
NGC 4374, NGC 4472 and NGC 4649. NGC 4261 also 
hosts an X-ray emitting AGN but it is highly absorbed, 
so that the thermal gas dominates the spectrum below 
^2 keV, enabling accurate gas te mperature determina- 



tion even at the smallest scales (§ 3.2; Zezas et al.||2005 
EI06). We excluded NGC 4374 since it has a very dis 
turbed X -ray morphology, in compa rison to the other 



systems ( Finogucnov &: Jones| 2001), making acc urate 
spectral dep r oiect ion very challenging. Although Chu 
t:azov et all (^008|) argued that the gas in NGC 1399 is 



relatively close to hydrostatic, we also excluded that sys- 
t em on account of its disturbed morphology at all scales 
(gcharf et al. 2005) and the presence of multi-phase gas 



( Puotc||2002rpumphrey fc Buotc| |2006|) , both of which 
may introduce systematic uncertainties into our mea- 
surements. We have already published the results for 
NGC 4649 (H08); this paper focuses on the three remain- 
ing objects (although we also discuss updated results for 
NGC 4649, having applied the refined analysis methods 
outlined in the present work). We list the basic proper- 
ties of these three galaxies in Table El 

We have previously published an X^-ray mass analy- 
sis for both NGC 4261 and NGC 4472 ( p06| ), and found 
that the gas appeared close to hydrostatic. However, our 
previous analysis focused primarily on the properties of 
the dark matter. Moreover we adopted a simple, ad hoc 
form for the stellar mass profile, which introduced con- 
siderable systematic uncertainty into the stellar mass-to- 
light (M/L) ratio measurement and clearly precludes the 
accurate determination of a central black hole mass. In 
this paper, we therefore present a new analysis for these 
galaxies, using an updated treatment of the stellar mass 
model, and the improved mass-determination methods 
outlined in H08. NGC 1332 is an X-ray luminous lentic- 
ular galaxy with a relaxed X-ray mor phologv and a cen- 
trally peaked gas temperatu re profile (Buote & Canizares 



1996| ; Humphrey ct al. 2004). Of the three galaxies, only 
NGC 1332 does not have a published Mbh measurement 
derived from dynamical methods. Both NGC 4261 and 
NGC 4472 have Mbh~ 5 x lO^M^, determined from, re- 
spectively, gas kinematics ( Ferrarese et al. 1996 ) and stel- 
lar dynamics (Gebhardt et al. 200£j^ 



3. X-RAY DATA ANALYSIS 

3.1. Data reduction 

The region of sky containing each galaxy was imaged 
by the Chandra ACIS instrument in the ACIS-S config- 



uration. The details of each observation are given in 
Table |[ The archival d ata were reprocessed using the 
CIAO 4.0 (beta 3) and Heasofi 5.3.1 software suites, 
i n conj unction with the Chandra calibration database 
(Caldh) version 3.5.0. To ensure up-to-date calibration, 
all data were reprocessed from the "level 1" events files, 
following the standard Chandra data-reduction threads^. 
We applied the standard correction to take account of 
the time-dependent gain-drift a nd char ge transfer inef- 
ficiency, as implemented in the CIAC\ tools. To iden- 
tify periods of enhanced background (which can seri- 
ously degrade the signal-to-noise, S/N) we accumulated 
background lightcurves for each dataset from low surface- 
brightness regions of the active chips, excluding obvious 
point-sources. Such periods of background "fiaring" were 
identified by eye and excised. The final exposure times 
are fisted in Table 0. For NGC 4261, we processed both 
datasets independently, and then merged the final prod- 
ucts (images, spectra, spe ctral response files, etc.) using 
the procedure outlined in H08 , which involves correcting 
for relative astrometric errors by comparing positions of 
the detected point-sources. 

For each galaxy we generated a full resolution image in 
the 0.5-7.0 keV energy-band and a corresponding expo- 
sure map computed at anenergy of 1.7 keV. Point sources 
were detected with the CI AC wavdetect task, which was 
set to search for structure at scales of 1, 2, 4, 8 and 16 
pixels, and supplied with the exposure-map to minimize 
spurious detections at the image boundaries. The de- 
tection threshold was set to 10~^, corresponding to < 1 
spurious source detections per chip. All detected sources 
were confirmed by visual inspection, and, for each, appro- 
priate elliptical regions containing approximately 99% of 
its photons were generated. In Fig |^ we show smoothed, 
flat-fielded Chandra images of each galaxy, having re- 
moved the point-sources from the images. To produce 
these images we first replaced all the photons within the 
source detection ellipse of each so urce with artificia l data, 
using the algorithm described in Fang et al. (2009). The 
data were then fiat-fielded with the exposure map and 
smoothed with a Gaussian kernel. Since the S/N of the 
image varies strongly with distance from the centre of 
the galaxy, the width of the Gaussian kernel was varied 
with distance according to an arbitrary power law. For 
each galaxy, the X-ray image shows a relatively relaxed, 
centrally peaked morphology within the central ^10 kpc 
(where the stars are expected to be the dominant mass 
component) . 

In Fi g we overlay the r adio contours from the NVSS 
survey (Condon et al. 1998| ) on the NGC 4261 image a nd 
those from the VLA FIRST survey ( [Becker et aI||l995D on 
the image of NGC 4472. In both galaxies the radio lobes 
appear coincident with slight depress ions in the X-ray 



brightness, which m ay be "cavities" (Biller et al. 2004 



Croston ct al. 2005 ) . Although we have previously shown 
( H06| ) that the gas density and temperature can be fitted 



with physically meaningful hydrostatic solutions even if 
data from the vicinity of these cavities is included, to 
minimize potential systematic uncertainties in our sub- 
sequent analysis we opted to mask them out, as shown 
in Fig |l|. We also masked o ut data around the jet in 
NGC 4261 ( [Zezas et al.| [20051 ) . 



http: / / cxc. harvard, edu /ciao/ threads / index, html 
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TABLE 1 
The galaxy sample 



Galaxy 


Distance 


Lj 




Age 


<[Z/H]> 


ObsID 


Date 


Exposure 




(Mpc) 


(1O"L0) 


(km s-i) 


(Gyr) 








(ks) 


NGC 1332 


21.3 


0.86 


321±14 


4 1+8-8 


0.32±0.30 


4372 


Sep 19 2002 


41 


NGC 4261 


29.3 


1.4 


308±6 


15±1 


-0.03 ±0.10 


834 


May 6 2000 


34 














9569 


Feb 12 2008 


101 


NGC 4472 


15.1 


1.7 


294±3 


9±2 


0.17 ±0.12 


321 


Jun 12 2000 


34 



TE. — Basic data for the objects in our sample. Distances were derived from the SBF survey of 



Tonr 



ot al l 



9nni )^xLu;'ecting the distance modukis by -0.16 mag to account for revisions to the Cephaid zero-point (e.g. Jensen 
pt al. 20031). Total J-band luminosities (Lj) were derived from our own modelling (§ Idl). The central velocity 



dispersion of the stars (ct*) was. 



from our Lick index analysis (Humphrey . 
(ObsID) 



taken from HyperLEJJA 



Buotc 



2006 



H06) 



and mean mctallicity (<[Z/H]>) were derived 
The Chandra observation identification number 



given for each dataset, and wc hst the date the observation was made and the total exposure time 
following removal of background flare events. 



NGC 1332 





Fig. 1. — Smoothed, point-source subtracted Chandra X-ray images of the galaxy sample in the 0.5—7.0 keV band (see text). These 
images only cover the central part of the ACIS-S3 chip, where the count-rate is sufficiently high for the structure of the X-ray emission 
to be clearly discerned. All three galaxies show relatively relaxed morphologies, although there is evidence of cavities in NGC 4472 and 
NGC 4261. The image contrast and smoothing scales were arbitrarily adjusted to bring out key features; the smoothing scale varies from 
~1" in the central part of the images to ~1' in the outer regions. We overlay radio contours from the NVSS survey on the image of 
NGC 4261 and VLA FIRST survey contours on the Chandra image of NGC 4472. We also mark regions of the field of view which were 
excluded from our analysis (shown circumscribed by dotted white lines). 

instrumental response may be uncertain. We extracted 
products from all active chips, excluding S4 (which suf- 
fers from considerable "streaking" noise). Appropriate 
count-weighted spectral response matrices were gener- 
ated for each annulus using the standard CI AC tasks 
mkwarf and mkacisrmf . For each spectrum, we esti- 



3.2. Spectral analysis 

We extracted spectra in a series of concentric, contigu- 
ous annuli, placed at the X-ray centroid of each object. 
We determined the centroid iteratively by placing a 0.5' 
radiu s aper ture at the nominal galaxy position (obtained 



from NED) and computing the X-ray centroid within it. 



The aperture was moved to the newly computed cen 
troid, and the procedure repeated until the computed 
posi tion converged. — The final ocntroido wore all 



mated the background using the method outlined in H06 
Spectral-fitting was carried out in the energy-band 0.5 



l^ Y^y^ 7.0 kcV. us ing the XSPEC spectral- fitting package (Ar 



0.6 1' of the optical centre , estimated from our HST data 
(§ ^), having corrected for differences in the absolute as- 
t rometry by matching X- ray and optical point sources 



naud 



1996 ) . We have shown previously that fits to Pois- 



( [Hui pphrcy fc Buoto || 20 9gD 



son distributed data which minimize x can yield sig- 
nificantly biased results, even if we rebin t he data to 
more than the canonical '^20 counts per bin (Humphrey 



I'fJ c widths of the annuli were chosen so as to con- 



tain approximately the same number of background- 
subtracted photons in each (we relaxed this criterion for 
the innermost bins to provide as high spatial resolution as 
possible), and to ensure there were sufficient photons to 
perform useful spectral-fitting. We placed a lower limit of 
2.5" on the annulus width (for the circular central bin, 
this refers to its diameter), to ensure the finite instru- 
mental spatial resolution does not lead to strong mixing 
between the spectra in adjacent annuli. The data in the 
vicinity of any detected point source were excluded, as 
were the data from the vicinity of chip gaps, where the 



et al J 120091) In contrast, we found that the C-statistic of 
pasf ( 1979 ) typically gives relatively unbiased results. 
We therefore performed the fits by minimizing C, but 
nonetheless took care to estimate any residual bias on 
our fits by using t he Mo nte Carlo technique outlined in 
Humphrey et al. (2009), as discussed in detail below. 
Although not strictly necessary for a fit using the C- 
statistic, we rebinned the data to ensure a S/N ratio of 
at least 3, and a minimum of 20 photons per bin. Such 
rebinning is useful during fitting as it emphasizes differ- 
ences between the data and the model if the fit is poor. 

For each galaxy, we fitted the data from all annuli si- 
multaneously. This allowed us to take advantage of the 
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NGC 4261 





Energy (ke7) 



Fig. 2. — X-ray spectrum measured in the central bin of each galaxy. We show the measured data (without removing instrumental 
features), and the best-fitting model, folded through the instrumental response (solid black line). We also show individually the modelled 
contributions of the hot gas (dashed red line), the unresolved LMXBs (dash-dot blue line) and, for NGC 4261, the central AGN (dotted 
magenta line). Given known correlations between the K-band light and the luminosity of unresolved point sources, we show the region in 
which the LMXB emission is expected to lie at l-cr uncertainty (upper, light blue shaded region), and the region in which the contribution 
from the composite X-ray emission of cataclysmic variables and stars is expected to lie (lower, grey, shaded region). Hot gas unambiguously 
dominates the spectra below ~2 keV in all three galaxies. 
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Fig. 3. — Density (p; triangles) and temperature (kT; asterisks) profiles for the hot gas in each galaxy. We also show the bias fraction fj, 
(see text) for both kT and p, which indicates that the derived data are largely unbiased. The best-fitting density and temperature models 
described in § M are shown as solid lines and fit the data well. We also show_(dotted lines) the best-fitting arbitrary parameterized models 
for the density and temperature used in our "traditional" mass analysis (§ 8.4). All error-bars correspond to 1-a uncertainties. 

pro.i ct model in Xsyec, which enables spherical depro- 



jecti( [)n. Spherical symmetry is typically assumed in X- 
ray mass analysis, even if the X-ray isophotes are not 
perfectly round, as deviations from this approximation 
are only expected to contribute a relatively small error 



into the recov ered mass profile (e.g. Piffarctti et al.||2003 
pavazzi 2005). We will return to the possible impact 
of asphericity in a future paper; for now, however, we 
are interested in how accurately a conventional spherical, 
hydrostatic mass analysis can recover the mass profile in 
the inner part of the galaxy. To take account of pro- 
jected emission from regions outside t he ou termost annu- 
lus we used the procedure outlined in H0(3| . To model the 
hot gas we adopted a vapec component, plus a 7.3 keV 
thermal bremsstrahlung component to account for un- 
detected low-mass X-ray binaries (LMXBs) (this model 
gives a good fit to the composite spectrum of the detected 



LMXBs in near by galaxies: Irwin et al. 2003 ; Humphrey 
fc Buotel ^0081 ). For NGC 4261, we included an addi- 
tional absorbed powerlaw compon ent in the central bin 
to account for the AGN emission ( HG6 ). We adopted a 
slightly modified form of the Xspec vapec code to en- 
able us to tie the ratio of each elemental abundance with 
respect to Fe between each annulus, while the Fe abun- 
dance (Zpo) was allowed to fit freely. To improve S/N, 
Zpc was tied between adjacent annuli where necessary. 
We allowed the global ratios of O, Mg and Si with respect 
to Fe to fit freely, while the remaining r atios were held at 
their Solar values ( |Asplund et"aIl|2GG4| ). For NGC 4472, 
we additionally freed S and Ni and for NGC 1332 and 
NGC 4261 we freed Ne. The absorb ing column density 
(Nh) was fixed at the Galactic value ( Dickey fc Lockman 
1990 ). The best-fi tting abundances were con siste nt wit h 
those reported in Humphrey & Buote ( ^006 ) and H06. 



6 



The X-ray spectrum of the innermost bin for each 
galaxy is shown in Fig ||, along with the best-fitting 
model, folded through the instrumental responses, and 
the contribution of each of the spectral components. 
Overlaid we show the expected spectrum of the unre- 
solved LMXBs (upper shaded region), and the expected 
integrated X-ray emission of cataclysmic variables and 
stars (lower shaded region). These were estimated from 
the relevant cor relations with the K - band light reported 



sev et al 



by, r espec t ively. Humphrey & Buote (2008) and Revnivt- 
( ^008 ) . Since neither correlation with the K- 



band light is exact, we shade the region allowed by the 
l-cr scatter in the relations. As is clear from the fig- 
ure, the measured normalization of the bremsstrahlung 
component is consistent (within ^1. 8-cr) with the LMXB 
flux p redicted from the relation of [Humphrey fc Buote 
(2008) for all three systems. The K-band magnitude was 
actually estimated from the I-band HST data (§ ^), and 
corrected for the difference in filters by using the differ- 
ence in magnitudes measured between the I-band data 
and K-band 2MASS data in a ~7" aperture. Clearly in 
all three galaxies, the X-ray emission in the innermost 
bin is completely dominated (below ~2 keV) by the hot 
gas emission. This is also true for the other annuli. 

Error-bars wer e com puted via the Monte-Carlo tech- 
nique outlined in HOG , and we carried out 100 error sim- 
ulations. We show the measured kT and gas density, pg 
(which can be computed from the normalization of the 
vapec comp onent) as a functi o n of r adius in Fig |^. As 
discussed in Humphrey et al. ( 2009 ), such simulations 
can be used to assess whether there is any bias on the 
best-fitting results by evaluating fb (the ratio of the bias 
to the statistic error) for each parameter. This is com- 
puted by taking the difference of the best-fitting value 
of each parameter and the mean value of that parame- 
ter obtained from the simulations, and then dividing this 
result by the statistical error. If |fb| -C 1 the fits are prac- 
tically unbiased. In Fig ^ we show fb for kT and pg in 
each annulus. We found that, in general, |fb| <C 1, sup- 
porting our use of the C-statistic. The temperature and 
density appeared slightly biased only in to two annuli in 
NGC4261, for which the error-bars are large and hence 
these data are unlikely to drive the temperature and den- 
sity fits discussed in § 0. The origin of the bias is unclear, 
but we suspect it may oe related to "deprojection noise" . 

4. THE STELLAR MASS DISTRIBUTION 

The modelling technique we have introduce d fo r 
analysing the X-ray data of early-type galaxies ( |H08| ) 
relies on an accurate (3-dimensional) representation of 
the stellar mass component. This is particularly true 
when one wishes to measure the mass of the central 
black hole, since one requires very accurate subtraction 
of the stellar mass. To determine the stellar light den- 
sity, we adopted an axisymmetric deprojection technique; 
this is particularly important given the highly ellipti- 
cal isophotes of NGC 1332, which is an almost edge-on 
lenticular galaxy. Although we subsequently spherically 
averaged the three-dimensional stellar density for use in 
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CVs and stars we use 



Revnivtsev et al 



(2008), for the spectrum of the 



abundances relative to Anders 

we hxed the ratio o 



A I. piggrri 



jstith kT=0.5, 



Solar 

Grcvesse| ( |l989| ), plus a E = 1.9 
the total flux in the 



our mass modelling (§ ^), we adopted this procedure to 
ensure the spherically averaged profile was as accurate 
as possible. Specifically we used the iter ative axisym- 



metric deprojection algorithm pioneered by Binney et al 



( |1990D 



This is a very general technique for axisymmetric 
deprojection, which is well-suited for modelling isophotal 
shapes ranging from highly elliptical (as in NGC 1332) to 
very boxy (as in the central part of NGC 4261), as well as 
systems in which the projected isophotal ellipicity varies 
stron glv as a function of rad i us, as is generallv the cas e 



here ( Kormendy et al. 2008; Buote fc Canizare^ 199(j ) 



powerlaw component 
2.0-7.0 keV band to the that in the 0.5-2.0 keV band to be 0.77 



As an axisymmetric method it cannot, however, account 
for the mild isophotal twists which are seen, so there will 
be some associated uncertainty in the deprojected light 
profile. Briefly, this technique involves first fitting an 
arbitrary function to the images, which is analytically 
deprojected onto an (r,0) grid as a "first guess" for the 
true stellar light distribution (here the light is assumed 
to be axisymmetrically distributed, with coordinate r in- 
dicating the distance from the centre of the galaxy and 
9 the zenith angle ) . A series of Lucy- Richardson "relax- 
ation" iterations ( Lucy 1974 ) subsequently improves the 
estimate nonparametrically. 

To obtain the optical light profile, we used 2MASS J- 
band images and, to provide superior spatial resolution 
in the crucial inner parts of the images, HST WFPC2 
data using the F814W filter (approximately I-band) for 
NGC 1332 and NGC 4472, and NICM0S2 data, using 
the F160W filter (approximately H-band) for NGC 4261. 
The NICMOS data were used in preference to WFPC2 
data for NGC 4261 to minimize attenuation due to the 
central dust disk, which complicates the deprojection. 
Since we aimed to model the temperature and data den- 
sity points measured with Chandra, for our purposes we 
required only the total stell ar li ght enclosed within a re- 
gion no smaller than ^1" (§ 3.2) and it was therefore not 
necessary to deconvolve the instrumental PSF from the 
HST data The prep a ration of the WFPC2 data is de- 
scribed in Humphrey ( ^009 ), although for NGC 4472 we 
used different data than those listed therein (specifically 
we used observation numbers u21g0603t and u21g0604t, 
for which the galaxy is better centred on the PC chip). 
For the NICMOS data, we use d the public l y avai lable, 
pre-cleaned image described in [Floyd et alj ( |200^ ). To 
prevent the emission from the weak AGN that is visi- 
ble in the NICMOS image of NGC 4261 from introducing 
bias, we replaced the data within a central ^0.35" region 
with the mean value of the surrounding pixels (which is 
reasonable given t he flat surface brightne ss profile in this 
part of the galaxy: Kormendy et al. 2008| ). In practice we 
did not use the entire WFPC2 and NICMOS images, but 
simply a portion close to the centre of the galaxy, chosen 
to overlap the area where seeing seriously compromises 
the 2MASS data. The images are shown in Fig |[ 

To obtain the initial guess for the stellar density, we 
simultaneously fitted models to the HST and 2MASS im- 
ages using custom-made software built around the MI- 
NUIT library^. We masked out obvious point-sources 
and the central parts of the 2MASS images, where 
the seeing is problematic. The WFPC2 images of 
NGC 1332 and NGC 4261 revealed central dusty discs; 

^ http:/ /Icgapp.cern.ch/project/cls/work- 

packages / mathlibs / minuit / index, html 
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Fig. 4.— Left panels: J-band 2MASS, I-band HST WFPC2 and H-band NICMOS images of each galaxy. For the HST data we annotate 
the image with the appropriate HST filter. Central panels: residual images, having subtracted our projected model for the three-dimensional 
light distribution (see text). Regions which were excluded from the deprojection are shown as white ellipses. Overall, these models provided 
a good fit to the optical light, except for the central part of NGC4261, where residual dust extinction is problematical. Right panels: the 
spherically averaged luminosity density as a function of geometric radius from our model (solid lines). For NGC4261 and NGC447?, wp. 



ly (dotted lines) rlpprnjectpH V-h 



1996h to the data of 



Kormcndy et al 



<Tid prnfilpg light profiles, arbitrarily scaled, obtained by applying the algorithm of 



Gebhardt et al 



(2008). These profiles agree well with those obtained using our default deprojection method. We 
also snow the best-fitting deprojected [Dersic model parameterization (dashed lines) for these profiles (see text). 

of changing the value of i in § ^. We aUowed the Ser- 
sic index, effective radius, scaUng factor a and position 



for NGC 1332, we masked this out as best we could. For 
NGC4472 and NGC 4261 we were able to fit the data 
satisfactorily with a model which is the analytical pro- 
j ection of a spherically sy mmetric Sersic density profile 
(Prugnicl & Simicn |1997 ), having broken the symme- 
try by arbitrarily scaling one of the coordinate axes; i.e. 
p*{x,y, z) = ap^^s{\/ x"^ + {ayY + z^). Here p*{x,y, z) is 
the luminosity density profile, a is an arbitrary scaling 
parameter and p*^^ is the deprojected Sersic density pro- 
file. In projection, this produces elliptical isophotes, the 
ellipticity of which depend on both the factor a and the 
inclination angle i. In general i is not known, and so 
we arbitrarily adopted i = 90°. We discuss the impact 



angle of the isophotes to fit freely, but constrained them 
to be the same for both the 2MASS and HST data of a 
given galaxy, taking into account orientation differences 
between the images. The centroids of the models were 
allowed to fit freely in each image, to allow for possi- 
ble errors in the absolute astrometry. We found there 
were two best-fit solutions which could not be distin- 
guished statistically — one with a > 1 (oblate) and one 
with a < 1 (prolate). By default we adopted the oblate 
model, but consider the effect of using the prolate model 
in § 0. For NGC 1332, the data required two such Sersic 
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models, plus an inclined exponential disc (which can be 
deprojected trivially, under the assumption of an arbi- 
trarily small but finite opening angle for the disc). We 
required the position angle of the Sersic models and the 
disc to be the same, and similarly tied the inclination of 
each component. Since the ellipticity of the disc model 
depe nds only on i. this model allowed us to constrain i 



73° 



to 

Having obtained the initial fits to the surface bright 
ness profile, and hence an initial guess for th e deprojected 



densit y, we used the relaxation algorithm of Binney et al , 
( |1990D to refine our estimates. Although the algorithm 
they discussed is relevant for a single image, it is straight- 
forward to generalize it to a stack of images with regions 
masked out. Since an image is simply a list of pixels, 
one always has the freedom to define a new "image" , ev- 
ery pixel of which maps uniquely to one of the pixels 
in the image stack. Not every pixel in the image stack 
need have a counterpart in the new "image" , so that one 
can readily mask out problem regions. Although this 
transformation complicates the definition of the projec- 
tion kernel, in principle exactly the same methods can 
be a pplied to deproject this new "image" . In practice. 



we a so subtracted the fitted background level from the 
count-rate in each pixel before mapping it onto the new 
"image" , resetting to zero any pixels which fall below 
counts, to prevent erroneously treating the background 
as a part of the source emission. 

After each relaxation iteration, the "improved" density 
profile was projected onto the 2MASS and HST image 
planes. We fitted this new "model" to the data, allowing 
the overall normalization, the scaling factor between the 
HST and 2MASS data and the background levels to fit 
freely. We continued iterating in this manner until the 
reduction in the fit statistic on any given iteration was 
less than 1, at which point we assumed convergence. In 
Fig ^ (middle column) we show the residuals (data mi- 
nus model) obtained from our deprojection technique for 
each galaxy. In general, the residuals are small, except 
for the central arcsec of NGC4261, where the dust disc 
is visible even in the NICMOS data. Due to the residual 
dust extinction, it seems likely that our deprojected pro- 
file for this galaxy slightly underestimates the true cen- 
tral stellar density; we i nves tigate the possible impact of 
this underestimate in § |7.2| . In the right hand column 
of Fig 1^ we show the spherically averaged density profile 
for each system, defined as the mean luminosity density 
within a thin, spherical shell of given radius. We also 
show (dashed lines) the best-fitting spherical Sersic ap- 
proximation to our deprojected profiles (the parameters 
of these fits are shown in the figure). These are shown 
to guide the eye since, by default, we use the Lucy de- 



projected profiles in our analysis (although in § 7.2 we 
investigate how significantly our results are changed if 
we use the Sersic profiles shown in Fig ^). For these fits, 
we find Ro=2.3 kpc, 3.0 kpc and 4.2 kpc, respectively for 
NGC 1332, NGC 4261 and NGC 4472, while the Sersic in- 
dex, n, is 4.2, 3.1 and 2.9, respectively. The values of the 
Sersic index we obtained for NGC 4261 and NG C 4472 
are lower than found by Kormendy et al. ( 2008| ), who 
reported values of 7.5 and 6.0, respectively, based on V- 
band surface brightness fits. The reason for the discrep- 
ancy is that we fitted the entire inner profile including 
the core, but did not extend the fit to large radii, while 



Kormendy et al. fitted out to very large radii, but ex- 
cluded the core region, which they considered to be an 
interesting departure from the model. 

To test further the sensitivity of our results to the de- 
projection method, for NGC 4261 and NGC 44 72 we also 
used the algorithm of Gebhardt et al. (1996) to d epro- 



ject the the V-ban d surface brightness profiles in Kor- 
mendy et al. ( ^008 ). We assumed that the stellar isoden- 



sity surfaces can be represented by similar, axisymmetric 
spheroids (with axis ratios 0.75 and 0.81 for NGC 4261 
and NGC 4472, respectively) and adopted an inclination 
angle of 90°. The resulting, spherically averaged den- 
sity profiles are shown (arbitrarily scaled, for clarity) in 
the right column of Fig ^ as dotted (blue) lines. Clearly 
there is good agreement with our best-fitting profiles over 
a wide radial range. The V-band profiles have a higher 
density outside '-~^10 kpc, which may arise from imper- 
fect sky subtraction in the 2MASS data. Since the stel- 
lar mass density only enters into our calculations through 
the total enclosed gravitating mass profile, which is dom- 
inated by the dark matter halo outside a few kpc ( H06| ), 
this does not pose a significant problem for o ur analysis. 
We di scuss the effect of using the deprojected Kormendy 



et al. profiles on our results in § 7.2 



5. MASS MODELLING 

The temperature and density profiles of gas in hy- 
drostatic equilibrium are uniquely determined by four 
factors — the mass profile, the profile of the entropy 

proxy (s = kTue where rie is the electron density) 
the gas pressure at a fiducial radius and the total gas 
mass enclosed within the fiducial radius. We choose suf- 
ficiently small a fiducial radius that we can assume the 
enclosed gas mass to be zero. One can, therefore, invert 
the problem and use the temperature and density pro- 
files to recover the mass. We have found the most robust 
way to achieve this is a "forward fitting" method, where 
physically motivated parameterized models for the mass 
distribution and the entropy profile are used to predict 
the density and tempera ture p rofiles, which are then di 



rectly fitted to the data ( H08 ; also see GastaldcUo et al 



2007b for some variations on this approach 

Assuming spherical symmetry, we modelled the en- 
closed mass profile as a stellar component (which was 
assumed to be proportional to the J-band lumin osity en- 
closed in a sphere of given radius) plus an NFW (Navarro 



t al. 1997) dark matter component, a central supermas- 
sive black hole, and the self-consistently calculated gas 
mass. We allowed the stellar mass-to-light (M/Lj) ra- 
tio to be fitted freely, along with the black hole mass 
(Mbh) and total mass and concentration of the dark mat- 
ter halo. We modelled the entropy as a broken powerlaw 
plus a constant, i.e.: 



s So sif{r) 



(1) 



where: 



\ lo'' kpc J \ 10 fcpc j > Sbrfc) (2) 

and the parameters sq, si, strk and /?2 were param- 
eters of the fit. Furthermore, we allowed the central gas 
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pressure to be determined by our fit. For a more detailed 



description of the modelling procedure, see |H08 . 
As is standard practice, we assumed that each tcmper- 



aturc and density data-point derived from our spectral- 
fits was a Gaussian-distributed, independent random 
variable. This enabled us to compare the models to the 
temperature and density profiles by computing the 
fit statistic (since the data were not Poisson distributed, 
wc do not expect the fi ts to be biased as described in 



Humphrey ct al. 2009). In H08 we computed the entropy 
profile from the temperature and density data-points, 
and fitted the models to the entropy and temperature 
profiles simultaneously. In the present work, however, 
we preferred to fit the models to the density and temper- 
ature as they are the quantities derived directly from the 
spectral fits. Clearly this is functionally the same thing; 
if a model fits the temperature and density profiles well, 
it will also fit the entropy. For NGC 4472, we only consid- 
ered data within ^ 3.2'(~14 kpc) as th ere is evidence of 
X-ray asymmetries outs ide this range (Irwin & Sarazin 
1996; Biller et al. 2004). Although we have previously 
shown that reasonable hydrostatic solutions can be con- 
struct ed ev en when including the Chandra data at larger 
radii (H06), we did not include those data here so as to 
minimize any potential source of systematic uncertainty. 
To obtain robust measurements of the fit parameters, 
rather than simply minimizing we adopted a Bayesian 
approach, as discussed below. 

5.1. The Bayesian prior 

As with all Bayesian analysis, the choice of prior is a 
matter of considerable import. While for some of the 
fit parameters (below) this choice is natural, for other 
parameters (the central gas pressure and the parameters 
characterizing the entropy profile), it is less clear. In 
those cases it is tempting to resort to a fiat prior, but this 
can be problematical. Since a prior which is defined to 
be fiat for a given parameter will no longer be flat under 
arbitrary transformation of that parameter (for example, 
taking the logarithm), even a "flat" prior adopted in ig- 
norance arbitrarily imposes our preconceived notions on 
the fit. Fortunately, this choice does not always have a 
significant impact on the derived results; in some cases, 
the data may actually place much tighter constraints on 
the parameter than any reasonable prior, and so a flat 
prior is as good a default choice as any. (Formally this 
is thp ca.se if the dvna,micaJ ra,nge of the pa,ra,meter af- 



lowei 1 by the likelihood function is sufficiently narrt 



that the Jacobean introduced by a reasonable parame- 
ter transformation is approximately constant). Further- 
more, if one is uninterested in the particular value of 
a given parameter (i.e. intending to marginalize over it), 
the choice of prior is also relatively unimportant provided 
that parameter has little or no covariance with any pa- 
rameter of interest. Nevertheless, it is unlikely that one 
will know that the problem satisfies these criteria before 
actually exploring the Bayesian posterior, and so a rea- 
sonable approach is to cycle through multiple priors to 
examine the impact of each choice on the results. Adher- 
ing to this principle, by default we adopted flat priors on 
these problematical parameters (i.e. those describing the 
entropy profile and the pressure at the fiducial radius), 
and subsequ ently investigated the effect of adopting al- 
ternatives (§ 7.1). 



For the other parameters the choice of prior is much 
clearer. A priori, we would expect that Mr h for each 
galaxy wou ld agree with the Mbh-c* relation ([Tremaine 
et al. 2002) and so we assumed a logi^Mbh prior which 



is Gaussian distributed about this relation, given the in- 
trinsic scatter in this relation (taken to be 0.25 dex) and 
the error on cr, (Table |l|). We discuss th e se nsitivity of 
our results to the form of this relation in § |7.l| . Addition- 
ally, we have also explored the effect of using a flat prior 
for Mbh which, although poorly motivated, is useful for 
disentangling to what extent the measured constraints 
are due to the prior and due to the data (§ ||)^. 

Given the measured age and mean abundance of the 
stellar population in each galaxy (Table [T]) , stellar popu- 
lation synthesis models can be u sed to predict the stellar 
M/Lj and its uncertainty (§ [8.21 ). Although this provides 
a natural prior for M/Lj, it is of interest to compare 
the population synthesis predictions to those obtained 
from the flts. Therefore, to avoid circular reasoning, we 
instead used a flat prior for M/Lj, but conflrmed that 
consisten t re sults were obtained with the more physical 
choice (§ 7J). In practice M/Lj is sufficiently well con- 
strained by the data that this choice has little effect on 
the results. 

If the galaxies were chosen entirely at random (irre- 
spective of morphological type or central velocity disper- 
sion), one might simply adopt the halo mass function 
predicted by dark matter simulations as the prior on the 
virial mass of the dark matter halo. However this is 
clearly not appropriate here; a better choice would be 
to restrict the allowed range of the total baryon fraction 
within the virial radius, based on the predictions of a 
theoretical model, as was done in H06| . Unfortunately, 
the physics of feedback is highly uncertain, making the 
results then sensitive to the assumptions implicit in the 
adopted model. Alternatively, one could employ an em- 
pirical rel ation linking the mass of the galaxy to its host 
halo (e.g. Conroy fc Wcchslei 2009). We investigate this 
option in § [7] but, in the interest of more generality in our 
results, by aefault we used a uniform prior for the loga- 
rithm of the halo mass in the range logioMdm — 10^^- 
lO^^'Mo. For a prior on the concentration, we adopted 
the theoretical concentration-halo mass relation for re- 



laxed halos found by Maccio et al. (2008). We used their 
relation for the "WMAPl" ACDJvl cosmology since, of 
the cosmologies they considered, that is most con sistent 
with observ ations of galaxies, groups and clusters (Buote 



ct al.||2007|) 



5.2. Exploration of the Bayesian posterior 

The exploration of the Bayesian posterior for a multi- 
parameter model is typically carried out with Monte 

* We note that another possible choice here, that of a prior which 
is fiat on log Mbh (which is the same as a Prob(MBH)=l/MBH 
in linear space), would be a poor choice since it would imply a 
priori that Mbjj is as likely to be found in, for example, the range 
lO^-fO^M© as the range lO^^-lO^M©. Since the data for two of 
the galaxies have little discriminating power at very low masses 
(FigTa), the use of this prior would unrealistically skew the black 
hole masses low and, crucially, will be sensitive to the lower mass 
cutoff. Conversely, a prior which is flat in Mbh produces a pos- 
terior probability distribution that is closer to that produced with 
the (more realistic) Mbh-""* prior. Since the data tightly con- 
strains the upper limit on Mbhi the results using the fiat prior are 
relatively insensitive to the upper mass cutoff. 
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Fig. 5. — Central part of the temperature and density profiles 
for the gas in NGC4472, illustrating the effect of adding a black 
hole. We show the best-fitting model with a non-zero Mbh (solid 
lines) and the best-fitting model with Mbh=0. The latter model 
is a poorer fit to both profiles in the central bin; the difference in 
between the two models is 5.5. Error-bars correspond to l-cr 
uncertainties. 

Carlo methods; we emplo yed the publicly a vailable^ 
MuhiNe st co de vers. 1.0 ( |Feroz et aT| ^008| ; |Feroz & 
Hobfon 2GG§| ), which implements a robust and efficient 
nested sampling algorithm (a Monte Ca rlo integratio n 
technique optimized for these purposes: Bkilling 2004). 
We have verified that we obtained consistent results with 
a Metropolis algorithm based Markov Chain Monte Carlo 
code. Since nested sampling integrates over a finite vol- 
ume, in the absence of hard boundaries for any param- 
eters (such as the restricted range of dark matter virial 
masses we considered) we restricted each parameter to a 
±10-tT region about the "best-fitting" value. This was 
estimated by first minimi zing with the Levenberg 
Marquardt algorithm (e.g. 
vides for each parameter 



Press et al 



30th the 



1992 ), which pro- 
best-fitting" value 



and an estimate of the statistical error, a. On com- 
pletion of the nested sampling run, we checked for self- 
consistency between the initial estimates of a and the fi- 
nal error-bars. If the latter were significantly larger than 
the initial guess, we re-ran the nested sampling algorithm 
with our refined error estimate, thus ensuring at least a 
~ ±10-(T region of integration for each parameter. 

On completion, the nested sampling algorithm returns 
a set of data-points which sample the posterior reason- 
ably well. One can use these data to construct a proba- 
bility density histogram for each parameter, marginaliz- 
ing over the other parameters. From this histogram, we 
found both the most probable values and a 90% confi- 
dence region, defined by stepping out along lines of con- 
stant probability density until 90% of the total proba- 
bility was enclosed. We found that the results do not 
depend sensitively on the width of the histogram used. 

6. RESULTS 

The model described in § || was able to fit the density 
and temperature profiles well and the best-fitting mod- 
els are shown in Fig ^. Strikingly, we do not see a sharp 
central temperature spike arising from the gravitational 
influence of the black hole, similar to that previously re- 
ported for NGC4649. There is a modest temperature 
rise in the centre of NGC 1332, but this is mostly a con- 
sequence of the peaked stellar density profile and is a 

® http:/ /www. mrao.cam.ac.uk/softwarc/multinest/ 



feature of our models even if Mbh= 0. Of all the galax- 
ies, only our model fits to NGC 4472 show any evidence of 
a higher central temperature in the presence of the black 
hole, and even then the effect is very weak, as shown in 
Fig a (which compares the best-fitting temperature and 
density models with and without a central black hole). 
We discuss the reasons for this difference from NGC 4649 



B.l 



The interesting marginalized parameters from our fit 
are given in Table ^j. For each object, the dark matter 
halo mass and concentration were degenerate with each 
other and we could not place interesting constraints on 
the virial mass (at least for NGC 1332 and NGC 4472) 
without applying an additional prior, suc h as a restric- 
tion on the allowed baryon fraction (see H0£). Rather 
than imposing such an ad hoc constraint, since neither 
of the dark matter halo parameters strongly correlate 
with Mbh or M/Lj, we simply marginalized over them 
and do not report them here. For NGC 1332, we were 
able to fit the data using a steep (/3i — 0.96) entropy 
profile without any evidence of a break; we could not 
constrain Sbrk and so we fixed it to a large value, out- 
side the field of view. For NGC 4261 and NGC 4472 the 
entropy profiles were similarly steep in the centre of the 
galaxy, but fiattened outside ~5-10 kpc. We show the 
entropy data-points, derived from the temperature and 
density data-points, plus the best-fitting entropy models 
in Fig p[ This characteristic shape for the entropy profile 
is similar to that fou nd in recent observations of X-ray 
brigh t galaxy groups ( Jctha et al. ^ 2007|: Finoguenov et al 



20G7| ; pastaldello et al.||20G7a| r pun et aL||2008| ), hutsjrk 
is at a comparatively smaller scale. 

The measured radial distribution of total gravitating 
mass is shown in Fig ^ for each galaxy, along with the 
contribution from each of the separate mass components. 
We found that the stellar mass component dominates 
within ~5-lG kpc, with a significant contribution of dark 
matter required to explain our observations at larger 
radii in NGC 1332 and NGC 4261 (for NGC 4472, the 
data we fitted did not extend to large enough radii for 
the influence of the dark matter to be easily discerned). 
These results are consistent with our previous analysis 
for NGC 4261 (and NGC 4472), despite the much sim- 
pler stellar mass profiles used in that work (HG6). 

Since the stellar mass is non-negligible in the inner- 
most bin, the constraints on Mbh for these galaxies are 
less tight than we were able to obtain for NGC 4649. 
Unfortunately, this makes the results depend, to some 
degree, on the choice of prior for Mbh- In Fig we 
compare the marginalized probability density functions 
for Mbh for the case where the Mbh-c* is used as a 
prior and for the flat prior case with the probability den- 
sity function of the Mbh-ct* prior itself. For NGC 4472, 
there is a clear peak in the probability density when the 
fiat prior is applied, allowing Mbh to be constrained to 
Mbh=(1-041J];^2) X iqSMq (90% confidence), and exclud- 
ing Mbh=0 at 98% probability (or "2.3cr")i°. Applying 
the Mbh-c* relation as a prior causes the probability 

This result is also obtained with a simple frequentist analysis: 
the difference in between the best-fitting model and the best- 
fitting model with Mbh=0 is 5.5 (which corresponds to a 98% 
confidence region). On the basis of an f-test, the significance of the 
black hole detection is therefore ~94%. 
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TABLE 2 
Fit results 



Galaxy 


xVdof 


Mbh 
(IO^Mq) 


M/L, 
(MqL-I) 


so 

(keV cm^) 


Si 

(keV cm2) 


/3i 


^brk 

(kpc) 


/32 


NGC 1332 
NGC 4261 
NGC 4472 


11.9/8 

22.9/12 

22.1/14 


'-'■^^-0.28 

(0.49) 
44+0-28 

(0.48) 
(0.69) 


1 16+''-l2 
J^-^O-0.14 

(1.15) 
1 621+°-''''^ 

(1.627) 

1 /I7Q+0.085 

l.^(Ci_Q 

(1.486) 


55+0-31 
(0.39) 

(0.558) 

1 .1+0.25 
^^••^-0.37 

(1.23) 


33.3tl^ 
(34.4) 

(78.6) 

66.ll"4 
(66.7) 


(0.951) 

-1 or,c+0.066 

(1.332) 

^■•^^-0.13 

(1.3) 


i3.4t^:^ 

^(^13) 
(5.32) 


47+o'2i 
'J-*' -0.31 

(0.54) 

-0.16 

(0.63) 



Nj:)TE. — The fit results for each galaxy, adopting the Mbh-'^'* prior for Mbh- We show x /dof for the best-fit model 
(§ m) and, for each interesting parameter, the most probable value and 90% confidence regions, marginalized over the 
other parameters. Since the set of marginalized values is not necessarily the set of parameter values which maximizes the 
posterior, we also list in parentheses the parameter values which correspond to such a "best-fit" . 



density peak to shift and narrow somewhat, but our con- 
clusions are not qualitatively changed. For NGC 1332 
and NGC 4261 the prior has a more significant effect on 
the results, however; below 3 x IO^Mq, the (Mbh-c*) 
prior completely dominates the shape of the probability 
density function, while the X-ray data largely determine 
its shape above ~4xlO^Mo. Put another way, the X- 
ray data alone provide a firm upper limit for Mbh for 
these two galaxies, while the Mbh-c* relation eliminates 
the portion of parameter space in which Mbh— *0, which 
is inconsistent with the apparent ubiquity of SMBHs in 
giant early- type galaxies. 

In Fig p we compare our constraints on Mbh for each 
galaxy with measurements using other methods. For 
NGC 44 72, we compare to the st ellar dynamical measure- 
ment of Gebhardt et al. ( 2009 , omitting a dark matter 
halo), for N GC 4261 wc use the gas disc dynamics mea- 



surement of Ferrarcse et al. (19961) and for NGC 1332 
for which there is no published dynamical Mbh mea- 
surement, we simply compare to Mbh estimated from 
the M-(T relation. We also include the results for a re- 
analysis of NGC 4649, applying the same Bayesian meth- 
ods and stellar light modelling used in the present paper 
(and the stellar dynamical analysis of Shen & Gebhardt 



p009| ) 



For completeness, we show the results using both 
choices of prior on Mbh- As expected, for NGC 4472 
and NGC 4649, the constraints are relatively insensitive 
to this choice. There is clearly agreement between the 
X-ray and stellar dynamical mass measurements. Partic- 
ularly striking is that, for NGC 4649, the error-bars ob- 
tained from the X-ray method are roughly half the size of 
those found by the stellar dynamical modelling. The X- 
ray measurement for NGC 4472 is marginally higher than 
the dynamically-determined Mbh. However, in the re- 



vised analysis of NGC 4649 by ^hen fc Gebhardt] (|2009|), 
Mbh increased b y a fa ctor ~2 from that obtained by 



P^bhardt et al.| ( |2003D . Whether this increase is due 
to inclusion of a dark matter halo in the newer models 
models or improved dynamical modeling will require ad- 
ditional studies. Nevertheless, it is plausible that there 
could be a similar revision for NGC 4472 when the dark 
matter halo is included, which would bring it into even 
better agreement with the X-ray data. 

For NGC 4261, the data unambiguously rule out that 
Mbh determined from X-rays substantially exceeds the 
gas disk dynamical measurement. Although the lower 
limit from our analysis is largely determined by the prior 
(so we cannot place stringent limits on any bias due to 



nonthermal pressure support, although see § pTl| ), there 
is no evidence from our analysis of a discrepancy be- 
tween the gas disk dynamics measurement and our X- 
ray results. For NGC 1332, our results constitute the 
first direct measurement of the SMBH mass, and so in 
Fig ^ we compare it to the Mbh-c* relation rather than 
a measurement derived from stellar dynamics. We find 
consistency; the X-ray determined Mbh clearly does not 
lie significantly above this relation. 

7. SYSTEMATIC ERROR BUDGET 

In this section we address the sensitivity of our re- 
sults to various data-analysis choices which were made, 
including the choice of prior. In most cases it is difficult 
or impractical to express these assumptions through a 
single additional model parameter over which one might 
hope to marginalize, and so we adopt the pragmatic ap- 
proach of investigating how our results change if the as- 
sumptions are arbitrarily adjusted. We focused on those 
systematic effects likely to have the greatest impact on 
our conclusions, and list in Table H how the "most prob- 
able", marginalized Mbh and M/Lj measurements for 
each galaxy are affected by each test. These estimates 
constitute a likely upper limit on the magnitude of each 
systematic uncertainty and we stress they should not be 
added in quadrature with the statistical errors. We out- 
line below how each test was performed, so those readers 
uninterested in the technical details of our analysis may 
wish to proceed directly to § |[ 

7.1. Priors 

As outlined in § ^.l| , the choice of priors on several 
of the parameters are arbitrary. We therefore consid- 
ered the effect of replacing each arbitrary choice with 
another, reasonable prior. Specifically, for each parame- 
ter describing the entropy profile, and for the central gas 
pressure, we switched from a flat prior on the parameter 
to a fiat prior on its logarithm. The effect, summarized in 
Table ||as APrior (entropy), is typically smaller than the 
statistical errors. Similarly, we have experimented with 
replacing the fiat prior on M/Lj with a Gaussian prior, 
the mean and sigma of which is determined from the stel- 
lar population synthesis models; see § 8.2. The effect of 
this choice is summarized as APrior(M/L) in Table |[ To 
assess the importance of the Mvir prior, we experimented 
with instead fixing the M^nr to a canonical value , base d 



on the empirical relation of jConroy fc Wechslcr ( |200 



which links the virial mass of a halo to the stellar mass it 
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TABLE 3 
Systematic error budget 



Test 


NGC 1332 


NGC 4261 




NGC 4472 




Mbh 


M/Lj 


Mbh 


M/Lj 


Mbh 


M/Lj 




(IO^Mb) 


(MoL"M 


(IO^Mb) 


(Mo 




(IO^Mq) 


(MbL"M 


Most probable 


n c:n + 0.41 
"■■'^-0.28 


1 1fi+0'i2 


U.44_0.24 






1 47+0.09 

J^-*'-o.ii 


APrior (entropy) 


±0.04 


±0.02 


+U.U9 

— 0.04 




-0.01 


+0.17 


-0.02 


APrior(M/L) 


-0.04 


-0.03 


-0.13 




+0.07 


+0.05 


-0.02 


APrior(Mvir) 


+0.11 


-0.01 


+0.04 




+0.00 


+0.24 


-0.02 


APrior(Gultekin) 


-0.04 


-0.03 


-0.10 




+0.00 


+0.10 


+0.01 


APrior(*lat Mbjj) 


-0.27 


+0.00 


-0.16 




+0.01 


+0.35 


-0.02 


AStars(geometry) 






-0.09 




-0.19 


+0.08 
-0.01 


-0.11 


AStars (deproj ection) 






-0.02 




-0.14 


-0.22 


-0.24 


AStars(Sersic) 


-0.09 


-0.13 


-0.22 




-0.06 


-0.15 


-0.15 


ZAUlVi 


-0.09 


+0.07 


-0.07 




+0.10 


+0.07 


+0.09 


ABackground 


+0.09 


-0.14 


+0.02 




-0.06 


-0.16 


+0.08 


ANh 


-0.07 


+0.03 


-0.03 




+0.01 


±0.03 


+0.03 
-0.02 


AXRB 


-0.10 


+0.04 


±0.02 




-0.02 


-0.15 


±0.02 


APlasma Code 


-0.08 


+0.03 


-0.01 




-0.06 


+0.25 


+0.04 


ABandwidth 


-0.11 


±0.03 


-0.04 




-0.02 


+0.33 
-0.16 


±0.07 


AStatistic 


-0.06 


-0.01 


-0.06 




+0.00 


+0.01 


-0.05 


ACentroid 


+0.08 


-0.12 








+0.02 


+0.00 


ADistance 


±0.02 


±0.09 


-0.07 




+0.17 
-0.15 


-0.03 


±0.06 


ARotation 


-0.08 


+0.09 


-0.05 




+0.01 


+0.01 


-0.01 



Note. — Estimate of the likely impact of potential sources of systematic errors in our analysis on 
our measurements of Mbh and the stellar M/L.j. We give the marginalized value and 90% statistical 
errors for each galaxy (top line). For a series of tests (sec text), the change in the marginalized value 
caused by an arbitrary change of the underlying assumptions in our analysis is given. We stress that 
these systematic errors should not be added in quadrature with the statistical errors. In most (but not 
all) cases, the systematic errors are no larger than the statistical errors. 



hosts. We computed the stehar mass from the Lj values 
given in Table 13 and the best-fitting stellar M/Lj ratio 
for each galaxy. The effect of this choice is summarized as 
APrior(Mvir). Finally, we investigated the effect of our 
Mbh prior by adopting the revised and updat ed Mbh 



(7* relation for early-type galaxies reported by Giiltckin 
et al. ( |2009D as the prior. The results (APrior(Gultckin)) 
indicate that this has little effect on our conclusions. In 
addition, we also report the result of using the (poorly 
motivated) flat prior on Mbh (APrior (Flat Mbh)), as dis 



cussed in § 5T. In general, adopting the flat prior had a 
larger impact on the derived parameters than any of the 
other choices. 

7.2. Stellar model 



As shown in H06, careful modelling of the stellar light 
is necessary to obtain an accurate measurement of the 
mass-to-light ratio. Similarly, since we did not typi- 
cally resolve scales at which the black hole is completely 
dominating the gravitating mass, we also expect Mbh 
to be sensitive to this modelling. For NGC 4261 and 
NGC 4472, both of which have uncertain geometry, we 
experimented with using a prolate model instead of an 
oblate model, and changing the inclination from 90° to 
45°. In NGC 1332, the geometry is less uncertain since 
the presence of the disc allowed us constrain the inclina- 
tion and symmetry axis. The effect of these choices is 
listed as AStars (geometry) in Table 0. 

In § ^ we compared our deprojectea stellar density pro- 
flle s to those obtain e d by u sing the deprojection method 



of Gebhardt et al. (1996), flnding general consistency 
within the central ^10 kpc. Nevertheless, to quantify 
any systematic uncertainties associated with the depro- 
jection method, we also used the V-band proflles for 
NGC 4261 and NGC 4472. Given the different fllters, we 



corrected the M/L ratios obtained based on the colours 
measured in a ^30" aperture , takin g the V-band data 



from Bandage & Visvanathan (1978) and computing the 
J-band magnitudes from the 2MASS data discussed in 
§ 0. We found that use of this profile did not lead to sig- 
nificantly different conclusions (AStars(deprojection) in 
Table 0). For all three galaxies, we also investigated the 
effect of adopting a simple Sersic approximation for the 
stellar light, as shown in Fig |[ This is useful as residual 
dust extinction in the centre of NGC 4261 likely leads 
us to underestimate the stellar density, so that the true 
density profile lies intermediate between the preferred 
deprojected profile and the simple Sersic parameteriza- 
tion. The impact of this test is listed as AStars(sersic) 
in Table |[ All three tests had a noticeable effect on the 
best-fitting results, comparable to the measured statis- 
tical errors. Nevertheless, these sources of uncertainty 
should not qualitatively affect our conclusions. 

7.3. Dark matter halo 



Although we showed in H06 that the dark matter ha- 
los of a sample of early- type galaxies were consistent with 
NFW profiles obeying the concentration-virial mass re- 
lation predicted by numerical str ucture formation simu- 
lations (see also Buote et al, 2007 ) , in some recent stellar 
dynamical analyses it has been argued that a "c ored log- 
arithmic" p rofile may provide a better fit (e.g. Thomas 



;t al. 2007). Therefore, to test the extent to wnich the 
modelling of the dark halo affects our results, we experi- 
mented with using the cored logarithmic profile in place 
of NFW. We found this had a slight effect (summarized 
as ADM in Table comparable to the other sources 
of error. We will return to the question of the optimal 
shape of the dark matter halo in a future paper. 
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Fig. 6.— pjEntropy profiles for each galaxy, derived from the density and temperature data in Fig|3| We also show the best-fitting entropy 
models (§ H), which fit the data very well. Error-bars correspond to l-o" uncertainties. 
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Fig. 7. — Radial mass profiles for each galaxy, derived from our fits to the Chandra data. The solid (black) lines indicate the total 
enclosed mass, the dashed (red) lines indicate the stellar mass, the dotted (blue) lines were the dark matter contribution, the dash-dot 
(orange) lines are the black hole mass and the dash-dot-dot-dot (magenta) line is the gas mass contribution. The grey shaded regions 
indicate the l-tr errn r on the total mass distribution. Overlaid are a set of data-points (and l-cr errors) derived from a more "traditional" 
mass analysis (§ 8.4), which generally agree very well with the fitted models (we stress the models are not fitted to these data-points, but 
are derived separately) . The arrows indicate the spatial scale corresponding to 1". 



NGC 1332 




NGC 4261 





Fig. 8. — Marginalized probability density distributions for Mbh for each galaxy. We show, as a solid line, the probability density when 
the Mbh-o"* relation was used as a prior. The dashed (red) lines are the probability density function with a fiat prior and the dotted 
(blue) lines are the Mbh-""* relation prior itself; both have been arbitrarily scaled for clarity. Also for clarity, we smoothed the distribution 
functions with a fourth-order Savitzky-Golay filter, which spans ~20% of the Mbh range for each system. For NGC 1332 and NGC 4261, the 
shape of the probability density function below ~ 3 X IO^Mq is mostly determined by the prior, while above ~ 4 X 10*Mq, the constraints 
are mostly given by the X-ray data. The X-ray data alone provide only upper limits for these two objects. For NGC 4472, the shape of the 
probability density function is much less sensitive to the prior. 
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Fig. 9. — Left panel: the derived black hole mass (Mbh) from the X-ray data, assuming a flat prior for MbHi shown versus black 
hole determinations from stellar dynamics (for NGC4472 and NGC4649; stars), gas disc dynamics (NGC4261; circle) or the M-cr relation 
(NGC 1332; square). Error-bars correspond to 1-cr uncertainties, although for NGC 1332 and NGC4261 we show the 3-(t upper limit 
(marked as an arrow). The dotted line indicates y = x. Right panel: The same, but using the Mbh-'''* relation as a prior on Mbh- Note 
that the comparison of the NGG 1332 Mbh measurement using this prior with the Mbh-""* relation, as shown in the right panel, should 
be treated with caution. We find good agreement between the masses measured with the different techniques. 



7.4. Spectral modelling 

We investigated several potential sources of systematic 
tainty in our opootral analyoia. — Firot, wo conoid 



unce: 
ered 



the impact of our background treatment, b)^ adopt- 



Brighenti et"aL||2009| ). For N GC 1332 we adopted t he ma - 
jor axis rotation profile from Dressier & Sandage ( 1983|), 
for NGC 426 1 we used the minor axis profile from Bender 
et al.| ( |1994|) and for NGC 4472 w e use d the major axis 



ing the spectra derived from the standar d background 
template events files distributed with the CIAO Caldl, 



having renormalized them to match the count-rate m 
the 9-12 keV band. This choice had a measurable ef- 
fect on our results, comparable to the statistical un- 
certainties (ABackground in Table Isl). We also inves- 
tigated how our results are affected oy changing Nh by 
±25% (ANh), replacing the bremsstrahlung component 
with a r=1.5 powerlaw or adjusting its temperature by 
±25% (AXRB), using a MEKAL rather than an APEC 
plasma (APlasma Code), or changing the energy range 
over which the fitting was performed to 0.4-7.0 keV, 0.5- 
2.0 keV or 0.7-7.0 keV (ABandwidth). Two further tests 
we carried out were using the fit statistic rather than 
the C-statistic (AStatistic), which had only a very small 
effect, and moving the centroid of the extraction regions 
by up to ~0.5", consistent with the accuracy we esti- 
mate for our centroiding procedure (ACentroid). This 
latter test was not carried out for NGC 4261, since the 
X-ray detection of the AGN makes the centroiding very 
accurate. 

7.5. Other tests 

Two final issues we investigated were the uncertainty in 
the adopted distance to the galaxies, and the possibility 
of gas rotation. For the form er, we adjust e d the distance 
by the 1-cr errors reported in Tonry et al.| (| 20Gl| ) and the 
effect is summarized as ADistance in Table yTTo test the 
impact of the latter, we allowed the gas to rotate along 
with the stars and included an additional nonthermal 
pressure compone nt to the equation of hydrostatic equi- 
librium, following Fang et al. (2009). We consider this 
to be likely an upper limit on the gas rotation; in none 
of the galaxies do we find the X-ray isophotes are much 
flatter than the optical light, as might be expected if the 
gas is rapidly rotating (for a more detailed discussion, see 



profile of Banchez-Blazquez et al. ( |2007| ). Assuming that 
the gas rotates as a rigid body, we find (unsurprisingly) 
that rotation increases the amount of enclosed mass we 
measure, but the effect (ARotation) is not larger than 
the statistical errors on M/Lj or Mbh- 



8.1. 



8. DISCUSSION 
Black hole mass measurements 



Under the assumptions of hydrostatic equilibrium, 
spherical symmetry, a single-phase ISM and a stellar 
mass component which follows the optical light, we have 
constrained the masses of the SMBHs at the centres 
of three galaxies. Com bined with our previous mea- 
surement for NGC 4649 (H08), this demonstrates that 
hydrostatic X-ray methods are a practical, competi- 
tive means for SMBH mass determination. We found 
that the masses determined from the X-ray method 
are in good agreement with those found from other 
techniques — sp ecifically stellar kinernatics for NGC 4472 
and NGC 4649 (iGebhardt et~all|2GG9l |2003|). g as disc dy- 
namics for NGC 4261 (Fcrrarese et al 1996), and the 
Mbh-ct* relation for NGC 1332. This maKs NGC 4261, 
NGC 4472 and NGC 4649 three of only a smaU handful 
of galaxies which have Mbh determined by more than 
one method. 

The agreement between our results and the SMBH 
masses obtained from other methods (especially for 
NGC 4472 and NGC 4649, for which the constraints do 
not strongly depend on the prior), provides support for 
the assumptions in our analysis. Key among these as- 
sumptions is that of hydrostatic equilibrium. While one 
might expect deviations from this approximation to be 
most prevalent close to the SMBH, since the ejccta from 
accretion episodes are known to interact with the the ISM 
at large scales, the local freefall timescale is shortest in 
the central region of the galaxy so that hydrostatic equi- 
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librium will quickly be re-established if the gas is stirred 
up. Furthermore, the physics of the coupling between 
an AGN and the surrounding ISM remains very uncer- 
tain, particularly on such small scales. We would ex- 
pect deviations from hydrostatic equilibrium most likely 
to m anifest themselves in nonthermal pressure s upport 
(e.g. Churazov et al. 2008 ; Zappacosta et al. 2006 ), lead- 
ing the X-ray method to underestimate the true Mbh ■ As 
is clear from Fig ||, the X-ray measurement for NGC 4472 
is slightly higher than the measurement from stellar dy- 
namics, while those for NGC 4649 are very close, suggest- 
ing that any such nonthermal pressure is small. While, 
for NGC1332 and NGC 4261, there is also no evidence 
that the X-ray determined mass is systematically under- 
estimated (since the upper limits imposed by the X-ray 
data are consistent with the gas disk measurement for 
NGC 4261 and the Mbh-ct* relation for NGC 1332), the 
data do not allow us to rule this out, since the lower limit 
on Mbh is largely determined by the prior. To investigate 
this in more detail there is a clear need for higher-quality 
Chandra data, which should allow more stringent upper 
limits to be placed on Mbh- Since the Mbh-c* rela- 
tion is not sufficiently accurate for this purpose, this will 
also require an accurate Mbh determination from stel- 
lar dynamics for NGC 1332. Despite these concerns for 
these two systems, it is unlikely that any deviations from 
hydrostatic equilibrium are restricted solely to the cen- 
tral bin, and deviations at larger radius should produce 
residuals in our fit to the temperature and density profiles 
and a systematically under-estimated stellar M/Lj ratio. 
As we discuss in jj ^.3|, the good fits to these profiles at 



larger radii, combined with reasonable M/Lj ratios for 
each galaxy in fact suggest that hydrostatic equilibrium 
is a reasonable approximation. 

Another crucial assumption in our analysis is that the 
ISM is approximately single phase at any radius. Al- 
though this approximation is not exact for N GC 4472, as 



evinc ed by regions of cool gas in its centre (|Caon et al 



i 



OOP ), the good overall agreement between Mbh deter 
mined from the X-rays and optical methods strongly sug- 
gests that it is, nevertheless, fairly accurate (implying a 
small filling factor for the cool phase). The dusty cen- 
tral discs in NGC 1332 and NGC 4261 similarly reveal 
departures from a perfectly single phase ISM. While, as 
discussed above, we cannot rule out the possibility that 
the X-rays may underestimate the mass for these two 
systems, provided the filling factor of the cool phase is 
small (for example, if the gas is restricted to this disc- 
like geometry) , the density and pressure of the hot phase 
obtained from our X-ray measurements will be fairly ac- 
curate, and thus sufficiently accurate for our purposes. 

The agreement we found between the Mbh measure- 
ments not only provides strong support for our X-ray ap- 
proach but also crucial, independent verification that the 
dynamical Mbh measurement techniques are accurate. 
As discussed in § ^, the stellar dynamical measurements 
for NGC 4472 and NGC 4649 are marginally lower than 
Mbh obtained from the X-rays, but this could arise from 
the omissio n of the dark matter halo in the stellar mod- 
elfing (e.g. ^hen fc Gebhardt| [2009| ) . Notwithstanding, 
this comparison suggests that deviations from triaxiality, 
uncertainties in the inclination of the galaxy and patho- 
logical orbital structure cannot give rise to additional 
systematic errors larger than a factor ^2 in Mbh for ei- 



ther galaxy (unless they are compensated for by finely 
tuned systematic errors in the X-ray determined mass). 
The similar agreement for NGC 4261 also suggests that 
any systematic errors in the Mbh determination are 
not significantly larger than the current statistic al er- 
rors. This is interesting since Ferrarese et al (1996) used 
sparsely sampled HST FOS data, while higher-quality 
STIS data sug gest that non-Keplcrian motions may be 
non- negligible (Noel-Storr ct al. 2003). Since the prior 
dominates the probability density function at low masses 
the agreement we find partially refiects the fair agree- 
ment between the gas dis k measurement of NG C 4261 
and the Mbh-c* relation of Tremainc et al. (2002), which 
was derived in part from this data-point. Still, it is very 
unlikely that this one point strongly infiuenced either the 
shape or scatter of their fit, so, combined with the tight 
upper limit on the mass from the X-ray data — which is 
also consistent with this datum — the agreement is very 
encouraging. 

In NGC 4649, the primary eflfect of the SMBH on the 
gas was to introduce a sha rp central temperature peak , 
similar to that predicted by Brighenti fc Mathews ( 199S| ). 
However, none of the three galaxies in the present study 
shows such a pronounced effect. This difference in be- 
haviour can be understood in part as a consequence of 
the black hole mass. In NGC 4649, the SMBH dominates 
the mass profile for much of the region enclosed by the 
central bin. Technically, the stellar mass enclosed within 
this bin is larger, but what is more important is the mass 
evaluated not at the outer radius of the bi n, but at an 
appropriate intermediate position (see e.g. Gastaldello 
et al. 2007b). Using an appropriate radius, the black 



hole actually dominates the mass, contributing 77lg% 
of the total (see Fig ^), whereas in the other galaxies 
Mbh is almost an oraer of magnitude lower and the 
stellar mass largely dominates (formally, it contributes 
18 ± 10%, 35 ± 15%, and 65^22% of the mass in the 
central bin for NGC 1332, NGC 4261 and NGC 4472, re- 
spectively). Correspondingly, therefore, the SMBH has 
a smaller influence on the gas. Another key difference 
between NGC 4649 and the other galaxies is the shape 
of its entropy profile, which is very flat within ^500pc, 
in contrast to the approximately s oa R proflles at small 
scales shown in Fig H. Since a lower central entropy tends 
to reduce the temperature of the gas in the centre, any 
temperature peak is suppressed in NGC 4472, NGC 4261 
and NGC 1332. 

One possible concern with our analysis is the depen- 
dence of the best-fitting Mbh measurements on the Mbh- 
ct, relation, which we used as a prior. In particular the 
Mbh-c* relation is determined in large part from stellar 
dynamical measurements which do not include dark mat- 
ter halos. Since this omission is able to bias the measured 
Mbh, the exact shape or scatter in this relation could 
be in error, particularly at the critical high-mass end, 
which may have an impact on our measured parameters. 
However, if significant, we would expect this to revise 
upwards the mass expected from the prior, which would 
serve only to restrict further the parameter space allowed 
(Fig H) and shrink our error-bars. Moreover, we ob- 
tained qualitatively the same results when the fiat prior 
was used in place of the Mbh-ct* relation (albeit we then 
only obtained upper limits on Mbh for two of the four 
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galaxies) , giving us confidence that subtle revisions to the 
Mbh-c* relation will not change our conclusions. Still, 
for a precision comparison of Mbh determined by differ- 
ent techniques, and, in particular, constraining the Mbh- 
a.if relation itself from X-ray data alone, deeper Chandra 
observations will be useful as higher-quality data gener- 
ally result in reduced sensitivity to the priors. 

Although we have demonstrated that hydrostatic X- 
ray mass analysis is potentially a powerful means to de- 
termine Mbh in gas-rich early-type galaxies, such mea- 
surements are contingent upon spatially resolved X-ray 
spectroscopy at the smallest possible scales. Given the 
limitations of Chandra, it seems likely that this tech- 
nique will only be applicable to a handful more galax- 
ies in the near future (§ |^). However, with the spatial 
resolution promise d by future X-ray miss ions, in particu- 
lar Generation-X (Windhorst et al. 2006D , such measure- 
ments should become routine. A factor ^5 improvement 
in PSF would likely correspond to ~2 orders of magni- 
tude increase in surveyable volume and hence number of 
accessible systems. Further, with such spatial resolution, 
it will become possible routinely to resolve the sphere of 
influence of the SMBH, making the Mbh measurement 
no longer dependent primarily on one data bin. 

8.2. Stellar M/Lj ratios 

In addition to the Mbh measurements, our hydrostatic 
models also enabled us to constrain the stellar M/Lj ra- 
tios, under the assumption that the stellar mass follows 
light. In Fig|l^, we compare our measured values to those 
predicted by three sets of single-burst stellar population 
(SSP) models, given the stellar age and abundanc e listed 
m T able [l|. In particular, we used the models of [Maras- 
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which arises due to the improved modelling of the stellar 
light in the prese nt paper; the galaxy was assumed to be 
spherical in HOS j , whereas it is appr eciablv flattened out- 
side a few kpc ( [Peletier et al 1990[ ), which we take into 
account with the deprojection method outlined in § ^. 

Although there is considerable scatter in the synthetic 
M/Lj values predicted by the SSP models, we find good 
overall agreemen t between th e predictions of the models 
using a Kroupa or [Chabriei" IMF and the results of our 
fits. In contrast, the models using the Salpeter IMF dra- 
matically over-estimate M/Lj. While this general trend 
is consi stent with the conclusions of our previous work, 
in H06 our measured M/L ratios were ty pically s ignifi- 
cantly smaller even than those assuming a Kroupa IMF. 



hnd good overall agreement between the measured M/Lj ratios 
and the models using a Kroupa or Chabrier IMF, suggesting that 
the gas is close to hydrostatic. 

In particular, for NGC 4472, we measured M/Lk to un- 
derestimate the SSP value by ~40%. The origin of this 
discrepancy, which has vanished in our updated analysis, 
can be easily traced to the simpler form for t he stellar 
density profile used in ou r earlier work. In II06 , we mod- 
elled the stellar light as a Hernquis^ (1990) profile, which 
significantly differs from the density profile derived in 
§ H outside ^--4 kpc. This difference was exacerbated by 
our not choosing the parameters for the stellar model in 
a self-consistent manner, thus introducing considerable 
additional systematic error into our older M/Lk mea- 
surement (see § 7.2 of H06). 

While the agreement m the lower panel of Fig |l^ 
is striking, it is unlikely that real early-type galaxies 
are monolithic, single-burst stellar populations. Indeed, 
the presence of line-strength and col our gradients (e.g. 
Peletier et al. 199C; Trager et al. 2000) indicates that the 
metallicity may vary radially and at least some ellipti- 
cal galaxies host mu ltiple-aged stellar populations (e.g. 



http://www.dsg.port.ac.uk/~marastonc. At the metallicity 
of these galaxies, the differences between the models with a blue 
and red horizontal branch are negligible. 



Rembold et al.| P005|) . Fitting SSP models to the Lick 
indices of such a system in order to determine a mean 
age and metallicity may yield a biased result, although if 
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one population dominates the stellar light, the bias may 
be small. Such a bias is most problematical for the age, 
which is the most important factor influencing M/Lj for 
an old system. Older populations typically have larger 
M/Lj ratios and so, since the mean stellar ages we esti- 
mated for NGC 4649 and NGC 4261 are both > 13 Gyr, 
the "true" M/Lj for these systems cannot be higher than 
our SSP predictions, for a given IMF. While this could 
in principle allow the X-ray data to be reconciled with a 
Salpeter IMF, for these two systems this would require 
the stellar light to be largely dominated by stars aged 
less than ~8 Gyr, implying a very large bias on the mea- 
sured age. Given the good agreement for all the systems 
between the me asured M / L.t ratios a nd the SSP models 
computed for a Kroupa or Chabrier IMF, a simpler hy- 
pothesis is that the bias is actually small and the Salpeter 
IMF is not correct. 

If the stellar population properties do vary radially, 
this can cause a violation of the assumption that the 
stellar mass strictly follows the light. Nevertheless, in 
the very red photometric band we adopt here the ef- 
fect of varying the metallicity is very slight (e.g. signifi- 
cantly chan ging the metal licity of a ~10 Gyr population, 
assuming a 
Marast 

licity gradient is not a concern. Furthermore, if an age 
gradient is present which might significantly affect the 
M/L ratio, it is likely to be give rise to a substantial 
colour gradient. In practice, all four galaxies show only 
very sli ght evidence for any such gradient, at least in V-I 
colour ( Poulain| 1988; Poulain & Nieto 1994, see also the 
good agreement between the shapes of the deprojected 
V-band and J-band profiles shown in Fig ^ for NGC 4261 
and NGC 4472). Therefore, the assumption that stellar 
mass follows J-band light is likely to be reasonable. 

We note that simple models of galaxy formation sug- 
gest that the condensation of gas into stars should pro- 
duce a gravitational reaction on the dark matter halo, 
leading to a cuspier dark matter prof ile through a pro- 



1 



roupa 



|2001| IMF and the SSP models of 
t^|2005fT?OTn^JVH]=0.35 to -0.33 causes only 
■o cliange in the stellar M/Lj ratio), so a metal- 



cess tei7n£d_^^adi^ati£c ( Blumenthal et al 

1986; Gncdin ct al. 2004| ). Such models increase the to- 



8.3) 



an IMF rad- 
or a mixture 



tal gravitating mass in the centre of the galaxy substan- 
tially, thus lowering the stel lar m ass-to-light ratio one 
measures for a given object ( |II06| ) when the models are 
applied. Our data clearly leave little room for a substan- 
tial reduction in the measured M/Lj ratios if we require 
consistency with SSP models. This is in clear contradic- 
tion with the adiabatic contraction models as currently 
implemented, unless they operate in a fin ely tuned con 
spiracy with nonthermal pressure (£ 
ically different from those discussed 
of multiple-aged stellar populations in each galaxy. This 
result is not surprising since similar difficulties for the 
adiabatic contraction models have been reported for disk 
galaxies ( Dutton et al. 2007| ; [Gncdin ct al. 2007 ) and re- 
cent simulations have found baryonic condensation to a 
have a significantly smaller effect on the dark matter halo 
( lAbadi et aL|^009|) . 

Finally, it is interesting to compare our measured stel- 
lar mass to light ratios with those obtained from stel- 
lar dynamical methods. Unfortunately, most dynami- 
cal studies to date have involved significant simplifica- 
tions, such as omitting a dark matter component, as- 



suming spherical symmetry or limiting the generality of 
the orbital structure, which may prevent a fair compar- 
ison with our work. Of the those in corporating a dark 



matter halo, Kronawitter et al. (200C) constructed spher- 
ical dynamical models for NGC 4472, but restricted the 
phase space distribution by applying an expansion using 
a set of known basis functions. They found M/Lb ^8-10, 
correcting to our adopted distance, which corresponds 
to M/Lj~2. 4-3, a ssuming a B-J colour for this galaxy 
of 3.1 (from NED), which is larger than our best-fitting 
M/Lj. However, the blue band they adopted is more 
sensitive to any contamination by a young stellar pop- 
ulation, so that the stellar mass may not follow B-band 
light as faithfully as the J-band light. Furthermore, the 
simplifications in their modelling and the use of a loga- 
rithmic (rather than an NFW) dark matter halo poten- 
tial (which can affect the best-fitting M/L; § 7.3) make 
the implications unclear. More general axisymmetric, 
orbit-based models that incorporate a dark matter halo 
(also using the logarithi nic potential) hay e been con- 
structed for NGC 4649 by phen fc Gebhaldt| ( p009|) , who 
obtained M/hv = 8.0 ± O.9M0 (roughly -70%higher 
than our measurement, when correcting for the galaxy 
colour). The authors were not able, however, to explore 
fully a range of inclination angles, which could affect the 



recov ered M/Ly ratio (e.g. Gavazzi 2005 ; Krajnovic et al 
2005 ). Still, we disc uss the possible implications of this 
discrepancy in § |8.3| . 

8.3. Hydrostatic equilibrium 

Our results suggest that the hot ISM of these galax- 
ies must be close to hydrostatic. Strong, highly localized 
deviations from this approximation would likely manifest 
themselves as significant residuals in our fits to the den- 
sity and temperature profiles, which are not generally 
seen (Fig H). More global deviations would result in a 
systematically mis-estimated pressure and consequently 
a bias in the measured enclosed mass, which is not borne 
out by the good agreement between the measured stellar 
M/L J ratios an d th e predictions of SSP models for a rea- 



sonable IMF (§ 8.2). The agreement we find between our 



Mbh measurements and those from different techniques 
supports the hydrostatic approximation extending into 
the central parts of the galaxies (albeit less emphatically 
on account of the larger error-bars) . 

In practice the selection of galaxies which are not sig- 
nificantly morphologically disturbed in the X-ray, at least 
over the scales under scrutiny (excepting possible cavi- 
ties, which we excluded from our analysis), eliminated 
galaxies in which hydrostatic equilibrium is most likely 
perturbed. Still, deviations from hydrostatic equilibrium 
could persist which do not obviously disturb the image, 
due to bulk or turbulent gas flows, magnetic fields or 
cosmic rays. Each of these effects will likely contribute 
nonthermal pressure support, causing X-ra y mass tech 



niqvics to underestimate the true mass (e.g. Nagai et al 



200' 



iraT 



Zappacosta et al 2006: Churazov et alfp008| ; ^ang 
200H| ). If such nontherinal pressure were prevalent 



in these galaxies it must be compensated for by errors in 
the SSP models themselves, which must underestimate 
t he n iass by almost the same amount. As discussed in 
§ ^.2| , at least for two of the galaxies, NGC 4261 and 
NGC 4649, this cannot be achieved simply by allowing 
a complex star-formation history, and requires a modi- 
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fication of the IMF from the Kroupa) or Chabriei form. 
Adiabatic contraction, or a similar effect, would also in- 
crease the gravitating mass in the central part of the 
galaxy and thus could offset the effects of such putative 
nonthermal pressure without requiring a modification of 
the IMF. In either case, however, the magnitude of the 
effect would have to be, serendipitously, balanced with 
the nonthermal pressure. 

Taking at face value the M/L ratios in ferred from the 
PEGASE SSP models that use a |Kroupa| IMF (which, of 
all the tested models, agree best with the measurements) , 
one can place an upper limit on the average global frac- 
tion of the pressure which is nonthermal (a) by fitting a 
model of the form: 

M/Lj{observed) = (1 - a)M/Lj{SSP) (3) 

To take into account the errors on both x and v axes, 
we use d a p rocedure similar to that outlined in Pres^ 
et al| ( |1992| ). We obtained a marginally acceptable fit 
(x^/dof=:7.4/3) with a < 0.02 (3-cr confidence region), 
indicating that no more than a few percent of the pres- 
sure support could be due to deviations from hydrostatic 
equilibrium. More conservatively, we estimated the ef- 
fect of inherent errors in the models by cycling through 
the results for the three different SSP codes. While 
this does not account for the possibility of a complex 
star-formation history, the failure to account for that ef- 
fect would mean that a is over-estimated (at least for 
NGC4261 and NGC4649, which provide crucial high- 
M/Lj leverage). The largest allowed v alue of a was 
~0.18 (at 3-cr significance) for Maraston's models. 

Assuming the SSP models are correct, these con- 
straints on nonthermal pressure are particularly impor- 
tant since two of the galaxies, NGC 4261 and NGC 4472, 
host active AGNs and associated radio lobes that have 
carved out cavities in the ISM. The very small value of 
a dramatically illustrates that the mere existence of cav- 
ities does not immediately imply that the ISM is out of 
hydrostatic equilibrium in regions away from the cav- 
ities (which were excluded from our analysis). Such 
stringent limits on deviations from hydrostatic equilib- 
rium are consistent with the sim ilar ~10-20% limit re- 
ported by Churazov et al. ( p008 ) for two galaxies which 
are manifestly more disturbed than our sample. Theo- 
retical structure formation models which produce mor- 
phologically relaxed-looking objects similarly suggest no 
more than ■^^ 25% nonthermal pressure due to turbulenc e 



( [Tsai et alj |l994 |Evrard et al.||l996| ; |Nagai et al.| pOOTI ) 
While we would not advocate the ro utine use of objects as 



morphologically dist urbed as M 84 (pinoguenov fc Jone; , 
^OOlD or NGC 4636 ( |Jones et al!l|2002|) for mass analysis 



(for s ome of the resulting issues see Brighenti & Mathews 
1997), taken together with our work, these results illus 



trate that there is little evidence, either observational or 
theoretical, that mild morphological disturbances in the 
X-ray image translate into errors on the derived mass 
profile that are much larger than the other systematic er- 
rors we assessed in this paper. This point speaks to the 
recent criticisms of Diehl fc Statlcr ( 2007 ), who argued 
that the lack of a correlation between the ellipticities of 
the X-ray and optical isophotes for a galaxy sample indi- 
cated the ubiquity of strong deviations from hydrostatic 
equilibrium. In part the lack of such a correlation was 
driven by the inclusion of a significant number of objects 



with X-ray features carved by AGN ejecta. As discussed 
above, while these disturb the X-ray morphologies, they 
may not lead to substantial violations of h ydrostatic equi 



l ibrium . Furthermore, as pointed out by Brighenti ct al 



(2009), highly subsonic turbulence and gas rotation (far 
below the level at which the enclosed mass inferred from 
hydrostatic techniques is significantly biased) can sub- 
stantially affect the X-ray e lliptic ity at the small scales 
studied by Diehl & Statler ( ^007 ), thus seriously under- 
mining the premise of their study. 

Notwithstanding the small amount of nonthermal pres- 
sure implied by our results, larger, albeit still modest, 
discrepancies have been reported between the mass in 
ferred from stellar dynamics and hvdr ostatic X-rav meth 



ods at large scales in sorae sys tems ( Shen fc Gebhardt 
2009; Romanowsky et al. 2008| ). It is unclear to what 
extent the discrepancies arise due to deviations from hy- 
drostatic equilibrium in the studied systems and to what 
extent systematic errors in the stellar modelling play a 
role (especially those due to triaxiality and the uncertain 
inclination of the galaxies) . ^hen fc Gebhardt argue that 
the stellar dynamics in one system over-estimates the 
mass inferred from X-rays, allowing room for nonthermal 
p ressure b ut only if the IMF is closer to a Salpete r than 



Kroupa form. Conversely, Romanowsky et al, argue 



that the mass inferred from globular cluster dynamics 
in another system, NGC 1407, underestimates the X-ray 
measurement. This is more difficult to reconcile with 
deviations from hydrostatic equilibrium, unless the gas 
is globally outflowing within the central '^20 kpc while 
maintaining a hig h X-r ay luminosity and smooth-looking 
X-ray isophotes (II06). Nevertheless, these studies can 
be used to place, entirely independently of our discussion 
of the SSP models, strict upper limits on the systematic 
errors expected for general use of the X-ray method. This 
would imply the bias on Mbh is no larger than ±0.3 dex 
but, as discussed above, we expect the actual systematic 
errors due to the hydrostatic approximation to be far 
smaller in the present paper. Ultimately direct observa- 
tional constraints on the ubiquity of turbulent nonther- 
mal pressure in the ISM of early-type galaxies and galaxy 
clusters should become possible with the advent of high- 
resolution n on-dispersive X-ray sp ectroscopy, as enabled 



by ^sfro-ff (Takahashi et al. 200^ ) and the I nternational 



X-ray Observatory {IXO: |WMtcet al.||2009[ ) 



8.4. Traditional analysis 

In this paper, we derived constraints on the mass pro- 
files of a sample of early-type galaxies by a "forward fit- 
ting" approach, wherein we parameterized the mass and 
entropy profiles and used them to derive density and tem- 
perature profile models to fit to the data. In this section, 
we discuss an alternative, more traditional approach that 
has widely been used to model the mass profiles of galaxy 
clusters. We demonstrate that it gives consistent results 
with our preferred approach, albeit with larger system- 
atic uncertainties. This method involves first fitting pa- 
rameterized models to the density and temperature pro- 
files, and then inserting these models into the equation 
of hydrostatic equilibrium. 



M(< r) = r- 



kT 



dlnp dlnT 
GiiniH V dlnr dlnr 



(4) 



to infer the mass (e.g. Lewis et al. 2003). Here G is the 
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universal gravitational constant, /i is the mean molec- 
ular weight of the gas ('^0.62) and mn is the mass of 
hydrogen. 

The problem with the traditional approach is that the 
mass is primarily determined by the derivatives of arbi- 
trary functions fitted to the binned density and temper- 
ature data. In principle there exists an infinite family 
of parameterized models, each having an arbitrary num- 
ber of maxima and minima in any given data-bin, that, 
when averaged over the bin, equally well fit the data 
but which correspond to different mass profiles (many of 
which may be unphysical). The fine detail of the true 
density and temperature profiles can be disguised by the 
binning, and may not be produced by an ad hoc, smooth 
parameterized model. This problem is particularly seri- 
ous when the temperature and density data-points have 
large errors since the data may appear adequately fitted 
by a very smooth model that does not, for any set of 
its parameters, resemble the true profile for the inferred 
mass distribution. In such circumstances, the range of 
mass profiles statistically allowed may not even bracket 
the true mass distributio n. These systematic errors (e.g. 
[Gastaldello et al.| ^007h| ) were the prime moti vation for 
the forward-fitting method we introduced in EI06 and 
HOS (and used in the present work), since those meth- 
ods enforce a physical mass distribution. Furthermore, 
the models employed in the present paper only involve 
adopting an ad hoc parameterization for the entropy pro- 
file; hydrostatic equilibrium requires this to rise mono- 
tonically, thus precluding the erratically rising and falling 
shapes which can confound attempts to parameterize the 
binned temperature and density profiles. 

Despite these concerns, the traditional method can be 
of use if the temperature and density profiles are suf- 
ficiently close to a smooth, parameterized form to pro- 
vide a (frequentist) consistency check on our more formal 
(Bayesian) methods. We therefore selected appropriate, 
albeit arbitrary, smooth models to fit the density and 
temperature data (described in detail in the Appendix), 
and hence d erived mass "data - points" from Eqn |4| as 



described in Gastaldello et al. ( 2007b ). The choice of 



model was not only dictated by its ability to match the 
temperature or density data, but also to generate a mass 
model which broadly resembles the expected form. The 
best-fitting temperature and density models are shown 
in Fig 0, and the derived mass "data-points" are shown 
in FigM Clearly, there is excellent agreement between 
the forward-fitting and traditional profiles, giving us con- 
fidence that the mass models adopted in § g are not seri- 
ously in error. Notable deviations between the two meth- 
ods are only seen in NGC 4261, outside ^^15 kpc, which 
likely reflects the limitations of our ad hoc density and 
temperature parameterizations. 

To quantify the ability of the traditional method to 
place constraints on Mrh, we fitted the set of mass "data 
points" shown in Fig |7| with a model comprising the 
dark matter halo, the stellar mass and the black hole 
(the gas mass is negligible over the radial span of the 
data; Fig 0). We obtained a good fit for all three galax- 
ies (xVdof = 5.1/4, 8.1/7 and 13.9/8 for NGC 1332, 
NGC 4261 and NGC 4472, respectively), and measured 
Mbh = (1-00 ± 0.39) X lO^M©, (1.2 ± 0.2) x IO^Mq and 
(0.32 ± 0.23) X IO^Mq, respectively, while the stellar 



TABLE 4 
BoNDi Accretion Rates 



Galaxy 


rA 
(pc) 


logioMbondi 

([Mo yr-i]) 


NGC 1332 
NGC 4261 
NGC 4472 
NGC 4649 


4^+36 
24+15 

16511 


_1 04+0.66 
of;+0.56 

— ^•^0_Q 57 

r, O + 0.66 

^■■^-0.65 
_1 41+0.29 
^■^^-0.33 



Note. — Accretion radii and Bondi 
accretion rates for each galaxy, com- 
puted from our fitted density and tem- 
perature models (see text). Error-bars 
correspond to 90% confidence limits. 

M/Lj is 1.05 ± O.O8MqL0\ 1.57 ± 0.03MqLq1 and 
1.40 ± O.O7M0Lq^, respectively. Although these mea- 
surements are broadly similar to those obtained from the 
forward-fitting method, it is clear that the error bars do 
not uniformly overlap between the two approaches. This 
is because the small error-bars obtained from the tradi- 
tional method actually mask the much larger systematic 
uncertainties discussed above. To gain an insight into 
the magnitude of these effects, we cycled through each 
of the density and temperature parameterizations given 
in the Appendix, and recomputed the mass profiles, pro- 
vided these models were able to fit the data adequately 
well. This cycling had a dramatic effect on the mass pro- 
files, with some models even implying unphysical nega- 
tive masses for some parts of the radial range. Fitting 
the mass models directly to the derived mass data points 
for each parameterization, we found that Mbh varied 
from its "best-fitting" value by lO^M©, 1.2x10^Mq and 
0.7 X lO^M©, respectively, for NGC 1332, NGC 4261 and 
NGC 4472. Similarly, M/Lj varied from its "best-fitting" 
value by 0.16MoLq\ 1.6MqL~^ and 0.35MqL~\ re- 
spectively. These large changes indicate that systematic, 
rather than statistical, errors dominate the traditional 
measurements for these galaxies. 

8.5. Bondi Accretion Rates 

Given the reliability of our three-dimensional gas tem- 
perature and density models (as inferred from the ac- 
curacy of our mass fits), we were abl e to estimate the 
Bond i accretion rate for each galaxy (J Pi Matteo et al ' 



2003| ; |Pellegrini||2005| ; [Allen et al.||2006D without relying 
on the ad hoc extrapolations sometimes employed. For 
adiabatic, spherically symmetric accretion from i nfinity 
onto anSMBH, the accretion rate is given by ( Frank 



taLl|l992| ): 



M, 



bondi ' 



Pg_ 



= 0.0126 



M 



BH 



(lO^Me 



yr 



(5) 



where Cg is the adiabatic sound speed of the gas and s is 
the entropy proxy (in keV cm^), which is conserved by an 
adiabatic flow. If we assume that an adiabatic Bondi flow 
exists inwards of some "transition radius" , one can then 
compute the accretion rate from Mbh and the entropy 
measured at any point in the flow (and without needing 
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to know explicitly the transition radius). Clearly this is 
only a crude estimate of the accretion rate since the cal- 
culation ignores radiative losses in the flow and neglects 
properly accounting for the boundary conditions at the 
transition radius. Furthermore, our models assume that 
hydrostatic equilibrium holds at all radii so that if such a 
Bondi flow exists it could introduce systematic errors into 
our calculation of Mbh from the X-rays. Nevertheless, 
as discussed in § such errors do not appear large. 

Allen et al.' (EOM) assumed that the transition to 
a Bondi flow occurs at the "accretion radius" , {— 
2GMbh/cI{oo), where Cs(oo) is the s ound speed far from 
the region of influence of the SMBH: [Frank et al] |1992D , 
that is, the point at which the gas flow becomes super- 
sonic in a conventional Bondi flow. For only one of our 
galaxies can we actually resolve (for NGC 4649 it is 
~ 160pc ~2") and the entropy profil e is, indeed, very 
close to adiabatic at this scale (H08). Still, the best 



the jet power is, at most, comparable to the lowest val 
ues found by 



Allen et al 



fitting So for the other galaxies is non-zero, suggesting a 
similar flattening of the entropy profiles at the smallest 
(unresolved) scales. We therefore found it convenient to 
use the asymptotic (r 0) value of s to evaluate Eqn 0. 
We list the most probable, marginalized val ues and 90% 
confidence regions for Mbondi in Table ^. Allen et al 
( ^006| ) also evaluated Mbondi for NGC 4472 and their re- 



ported value agrees, withi n errors, wit h our result. 

For a sample of galaxies, Allen et al. (2006) found that 
Mbondi correlates tightly with an estimate of the total 
AGN power ejected in the jets. This jet power was esti- 
mated under the assumption that the thermal energy of 
the ISM is increased by the work done in creating cav- 
ities with a relativistic gas, so that the mean power in 
the jet can be obtained from the gas pressure, t he cavity 
size and an estimate of the time to create it. Mathews 



Brighenti (2008) argued, however, that cavities cre- 
ated with cosmic rays do not heat the gas in this way 
and should, instead drive large-scale circulation of the 
gas. In this case, the jet power calculation is likely a 
very inaccurate way to estimate the energy injected into 
the ISM by the AGN. Despite the se concerns, the tight 



correlation found by Allen et al. is intriguing, and it 



is interesting to investigate whether the galaxies in our 
sample are consistent with it. 

The values of Mbondi in Table ^ are comparable with 
the mediuin -to-high accretion rate systems reported by 



Allen et al., implying that the jet power should be con- 
siderable for all of these systems if they lie on the same 
relation. While there is clear evidence of AGN-blown 
cavities in NGC 4261 and NGC 4472, the situation is 
much less clear for the other two systems, which have 
higher Mbondi- The X-ray image of NGC 1332 is smooth 
and relaxed and exhibits little or no evidence of X-ray 
cavities (Fig IT]). Nonetheless, the modest X-ray surface 
brightness of this object, coupled with the lack of obvi- 
ous radio lobes (given the weak emission in the NVSS 
image) make the identification of small, low-contrast X- 
ray depressions challenging. More significant, however, is 
the absence of cavities in a d eep Chandra observation of 
the X-ray b right NGC 4649 (|H08| ) . Based on much shal- 
lower data, Shurkin et al.| ( 200§| ) claimed to find small 
X-ray cavities (which we did not confirm), but even tak- 
ing these at face value, the implied jet power estimate 
is only ~ 1.3 x lO^^erg s~^. For this system, therefore. 



while Mbondi compares to their 
highest Bondi rates, in stark discrepancy with their cor- 
relation. The existence of g alaxies with high Mbondi but 
no cavities suggests that the Allen et al. correlation takes 
time to be established; we may be observing NGC 4649 
and NGC 1332 at an unusual time, while the AGN is just 
turning on and the radio lobes are beginning to inflate. 

9. CONCLUSIONS 

We have presented new hydrostatic. X-ray models for 
the centres of three early-type galaxies observed with 
Chandra. Combined with o ur recent study of the ellipti- 
cal galaxy NGC 4649 (HO^), we found, in summary: 



1. The black hole masses measured by our method are 
in agreement with those obtained by stellar or gas 
dynamics techniques. This provides support not 
only for our new X-ray approach, but also for the 
accuracy of the dynamical methods. 

2. Like in NGC 4649, the black hole in NGC 4472 is 
unambiguously required by our mass models; even 
with a conservative, flat prior for Mbh, it is de- 
tected at the ~2.3-(T level. 

3. Accurate stellar stellar M/Lj ratios require careful 
modelling of the stellar mass component; with our 
detailed models, we obtained stellar M/Lj ratios in 
agreement with the predictions of single-burst stel- 
lar p opulation s ynthesis ( SSP) models computed 
for a Kroupa or Chabrier IMF. This leaves little 
room for the steepening of the central dark matter 
density profile predicted by models of "adiabatic 
contraction" , unless it exists in a conspiracy with 
nonthermal pressure. 

4. Taking the SSP models at face value, this agree- 
ment suggests the gas is very close to hydrostatic 
and nonthermal pressure can provide no more than 
^^10-20% of the total support, unless there is a con- 
spiracy between the shape of the IMF and nonther- 
mal pressure. 

5. The two galaxies with the highest Bondi accretion 
rates exhibit little o r no evidenc e of X- ray cavities, 
suggesting that the Allen et al. (2006) correlation 
with the AGN jet power takes time to be estab- 
lished. 
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TABLE 5 
Traditional analysis models 



Galaxy 



Model 



{g cm~^\ keV) 



M 

{g cm~^; keV) 



Rc 

(kpc) 



Rc2 
(kpc) 



(4.6 ±0.7) X 10"^=' 
(1.6 ±0.1) X 10-24 
(7.9 ±1) X 10-^^ 



Density models 



NGC 1332 Eqn 
NGC 4261 Eqn 
NGC 4472 Eqn 



4J 



(2.1 ±0.4) X IQ- 
(6.4 ±3.4) X 10- 



-27 
'27 



0.25 

0.198 ±0.008 
0.25 ± 0.02 



0.98 ±0.07 



42 ±8 
15± 18 



0.47 ±0.008 
0.543 ± 0.007 
0.48 ± 0.02 



NGC 1332 Eqn 
NGC 4261 Eqn 
NGC 4472 Eqn 



0.39 ±0.02 
0.609 ± 0.006 
0.65 ±0.01 



0.3 ±0.06 
0.75 ±0.06 
0.42 ± 0.05 



2.1 ± 1.4 
6.4 ±0.6 
4 ±0.5 



5 ± 0.5 
2.8 ± 0.4 
2±1 



Note. 



Best-fitting 



vaults for ttie arbitrary, parameterized models to the temperature and density profiles used in our 
"traditional" analysis (§ S.41. For each galaxy, we identify which of the functional forms given in the Appendix are used 



( "Model" ) and, for each rnoacl parameter, the best-fitting values and 1-cr errors. 
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APPENDIX 
TRADITIONAL ANALYSIS MODELS 

We here o utli ne the explicit models used to parameterize the density and temperature profiles in our "traditional 
analysis" (§ B.4). These models are entirely ad hoc and were motivated only insofar as they reasonably capture the 
global shape of the appropriate profile. In practice, the particular combination of density and temperature model 
preferred to fit the Chandra data of each galaxy (which is given in Table |) was chosen as that which most accurately 
reproduces the mass profile derived in § g. 

To fit the density profile we used either a "cusped beta model" , 



p = Ai 



or a "double beta model" , 




-3/3 



1 



-3/3 



and to fit the temperature data we adopted one of the following arbitrary models: 

R 



fcT = Ai + ^2 1 



Rr 



(Al) 



(A2) 



(A3) 



fcT = Al + As 1 



R 



Rc 



(A4) 



tic 



exp 



Rc \Rc 



The best-fitting parameters for each galaxy are given in Table |. 



1 — exp 



R 

R'c 



(A5) 
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